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ABSTRACT 

The  applicability  of  the  methods  of  system  identifi- 
cation to  the  determination  of  ship  model  motion  parameters 
is  studied  taking  as  an  example  a  model  of  the  Mariner 
class  hull  form. 

A  mathematical  expression  describing  the  dynamical 
behaviour  of  an  ocean  vehicle  is  first  presented.  Then  the 
values  of  the  ship  model  motion  parameters  to  be  identi- 
fied are  predicted  using  available  techniques.  Data 
representing  the  dynamical  behaviour  of  the  hull  form 
under  controlled  conditions  Is  obtained  with  model  tests 
in  the  towing  tank.  This  data  is  then  used  together  with 
the  first  estimated  values  of  the  motion  parameters  to 
determine  the  corrected  values  of  those  parameters. 

The  sensitivity  of  convergence  to  several  noise 
factors  caracteristic  of  the  process  is  also  studied. 

The  results  indicate  a  fairly  strong  dependence  on 
those  noise  factors  and  on  the  initial  uncertainty  that 
in  general  decide  the  final  convergence.  This  fact  seems 
to  indicate  that  the  particular  technique  used  and  the 
particular  input  adopted  are  not  in  conjunction  the  most 
appropriate  for  this  type  of  dynamical  behaviour. 
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CHAPTER  1 


INTRODUCTION 

Consider  a  system  whose  inputs  are  known  and  the  outputs 
observable.   System  identification,  in  general,  has  the 
objective  of  determining  the  structure  of  that  system 
using  the  known  inputs  and  measuring  the  outputs  in  order 
that  one  may  be  able  to  predict  the  behaviour  of  the 
system  (outputs)  to  the  various  excitations  expected  from 
the  environment  (inputs).   One  such  a  system  of  interest 
in  naval  architecture  is  the  ocean  vehicle.   The  naval 
architect  wants  to  be  able  to  predict  the  velocities, 
accelerations  and  other  important  motion  responses  that  a 
ship  will  experiment  at  sea  in  order  to  effectively  design 
por  proper  response,  structural  integrity,  and  other 
critical  aspects  of  operability . 

The  identification  of  a  system  of  completely  unknown 
structure  is,  in  most  cases,  impossible,  at  least  for  the 
time  being.   There  ere,  however,  some  types  of  systems  for 
which  identification  is  possible.   Those  consist  of  para- 
meter identification  problems  where  the  general  structure 
of  system  is  known  but  the  particular  values  for  a  set  of 
parameters  are  to  be  determined.   The  rc^^n    vehicle  has, 
fortunately,  these  characteristics.  'Using  rigid  body 
dynamics,  one  is  able  to  develop  the  so-called  equations 
of  motion  which  describe  the  dynamical  behaviour  of  the 
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vehicle.   Tn  these  equations  there  are  some  parameters 
(ship  motion  parameters  or  hydro-dynamic  coefficients)  which 
are  not  known  or  at  host  whose  value  is  known  with  a  certain 
degree  of  uncertainty.   It  will  be  necessary  then,  ir  order 
to  fully  predict  the  behaviour  of  ocean  vehicles  at  sea, 
to  obtain  the  correct  values  of  those  parameters  or,  at 
least,  by  selected  measurements,  to  decrease  the  degree  of 
uncertainty.   Unfortunately  these  measurements  are,  more 
often  than  not,  combined  with  noise  from  the  instruments 
themselves  or  from  the  environment  at  the  time  the  measure- 
ments take  place  or  from  both.   These  types  of  noise  need 
+  n  v.rs  taken  ^nto  account  when  identifvinf  the  Daram  ;i  ^>~~. 
An  examole  of  the  fo^efoine  is  the  c^se  of  experiments  in 
a  towing  tank  to  obtain  data  from  models  to  apply  to  full 
scale  vehicles.   Measurement  noise  exists  whenever  one 
uses  any  type  of  transducers  that  always  show  some 
hysteresis  and  are  characterized  by  a  specific  dynamical 
behaviour.   Environment  noise  will  be  caused  among  others 
by  interference  between  the  model  end  undesired  waves, 
and  vibration  of  test  equipment.   Of  course  these  effects 
could  be  included  in  the  mathematical  model  o+'  the  vehicle 
but  it  would  make  the  analysis  of  its  behaviour  much  more 
difficult.   Thus  it  is  seen  that  sometimes  it's  preferable 
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to  assume  a  simpler  model  for  the  behaviour  of  the  system 
and  introduce  the  effects  not  accounted  for  in  the  mathe- 
matical model  as  noise  (environmental  noise)  whose  general 
nature  is  assumed  known. 

Many  of  the  available  techniques  of  system  identifi- 
cation were  developed  for  a  dynamical  system  whose 
behaviour  may  be  expressed  in  a  so-called  state  space 
repres  entation^l^j : 

X   -   f(x,u,t)  ;    x(t0)  =      x0         (1.1) 
v  =   h(x,u,t)  (1.2.) 

where  the  vector  x  expresses  the  state  of  the  system  at 
time  t  (t>  tn)  ,  the  vector  y_  represents  the  output  of  the 
system  at  time  t  and  the  vector  u  the  incut  at  time  t. 
x   is  the  state  of  the  system  at  time  t0  (initial  conditions) 
and  x  is  the  derivative  of  x  with  respect  to  time'.   f  and 
hare  vectors  that  describe  the  structure  of  the  system. 
The  main  assumption  to  be  made  at  this  stage  about  the 
above  system  in  order  that  meaningful  results  can  be 
obtained  is  that  there  v'll  be  unique  values  of  x(t)  a.nd 
v(t)  for  a  given  u(t)  and  x(t0)  =  xQ  for  all  t  such  that 

The  purpose  of  any  technique  of  system  identification 
will  then  be,  by  operating  on  u(t)  and  y_(t),  the  deter- 
mination of  f  and  h.   As  said  before  it  will  be,  in 
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general,  not  possible  to  obtain  the  solution  to  this  pro- 
blem.  Therefore,  the  assumption  of  a  known  structure 
except  for  a  set  of  unknown  parameters  is  marie  and  the 
system  is  expressed  by  the  following: 

x  =  f(x,u,jD,t)         x(t0)  =  x0        (1.3) 
X  =  £(x»HfEt"t)  0  -h) 

where  the  vector  p_  represents  the  set  of  parameters  to 
be  determined . 

In  order  that  the  identification  of  those  parameters 
be  of  some  value,  it  is  further  assumed  that  p_  is  time 
invariant  in  the  time  interval  of  interest  t.^  t  <  t ,  or 
at  least  that  its  time  dependence  is  Vno'"r  meaning  that 
once  some  p_("^)  ,   t  <^<t-   ,  is  determii  ed,  p_(t)  will 
be  known  in  the  interval  [  t  ,  t£}. 

It  must  be  emphasized  at  this  stage  that  the  fore- 
going does  not  assume  the  uniqueness  of  that  representa- 
tion.  Only  the  uniqueness  of  the  input/output  relation- 
ship is  assumed.   That  is,  there  will  be  'for  a  given 
u(t  , t)  and  v(t),  a  unique  b  such  that 

v(t)  -  b(u(t0,t)  ,  t)  (1.5) 

It  will  be  possible  to  show  the  existence  of  a  repre- 
sentation x*  t    £' »  h1  end  p_'  such  that 

x»  -  f'(x',u,p_\t)         x-,(.t0)  -  XQ'    (1.6) 
v  =  h'ix'u.E'  ,t)  (1.7) 
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The  problem  is  then  redefined  as  the  determination 
of  a  set  of  parameters  p_  of  some  state  representation  of 
the  system  defined  by  v_(t)  -  b(u(tQ,t),  t). 

Coming  back  to  the  case  of  an  ocean  vehicle,  it  will 
be  shown  that  the  resulting  equations  of  motion  may  be • 
regarded  as  a  state  representation  of  the  vehicle 
dynamical  behaviour  end  that  the  assumptions  made  so  far 
are  met  by  that  mathematical  model. 

The  foregoing  assumes,  apart  from  the  unknown  set 
of  parameters  p_,  no  other  types  of  uncertainty.   But,  as 
said  before,  it  is  sometimes  advantageous  to  select  a 
simpler  math  5matical  model  for  a  given  system  and  intro- 
duce the  so-called  system  or  environment?!  noise  since  it 
results  from  interactions  between  the  system  and  the 
environment.   Also  consideration  must  be  given  to 
measurement  noise.   A  model  that  takes  into  account  those 
factors  will  be  of  the  following  form  [2]  : 

i   =  f(x.u.E.w.t)         x(t0)  -  xQ       (1.8) 
Z  ~-   h(x,u,p.,v,t)  0  .9) 

where  the  vector  w  represents  the  environmental .noise  of 
the  system  end  the  vector  v  the  measurement  no*?'0. 

Errors  in  the  structure  of  the  system  will  come  through 

as  environmental  noise. 
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Tn  addition  (case  of  time  varying  systems)  the  parameter 
vector  o  may  have  the  following  structure: 

E  -  g(E,t,w»)  o(t0)  =  e0     (1.10) 

where  g  is  known  and  w'  is  noise. 

In  order  to  simplify  the  above  formulation  the 
following  assumptions  are  generally  made: 

1.  The  mathematical  model  is  time ' invariant . 

2.  System  structural  uncertainty  will  appear  as 
an  apparent  additive  noise  process. 

3.  Measurement  noise  may  also  be  considered  an 
additive  noise  process. 

'4-,   Output  measurements  are  linear  functions  of 
the  state  of  the  system  and  are  structurally  independent 
of  the  input. 

5.   Parameters  are  coefficients  in  the  model 
structure  that  are  to  be  identified. 

The  system  dynamical  behaviour  may  be  expressed 
as  follows; 

x  -  f(x,u,p_)  -:-  w  (1.11) 

£  =  o  (1 .12) 

Y.  -  H  x  +  v  (1-13) 

where  H  is  the  measurement  matrix,  often  a  diagonal 

matrix. 
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The  purpose  of  the  id  entif  ioation  will  be  the  deter- 
mination of  p_. 

In  Chapter  2  the  equations  of  motion  for  an  ocean 
vehicle  will  be  developed  and  it  will  he  shown  that  indeed 
they  constitute  a  state  space  representation  of  the  ocean 
vehicle's  dynamical  behaviour.   It  will  be  shown  also  that 
they  satisfy  the  foregoing  assumptions. 

In  Chapter  3  "the  methods  of  predicting  the  ship 
model  motion  parameters  will  be  outlined  and  values  of 
those  parameters  for  the  mariner  class  hull  form  will  "be 
presented  (Appendices  1,  2,  3). 

In  Chapter  ^   the  excitation  from  the  environment 
w ill  be  d e s  c r i b ed  and  its  s  e v e r a 1  c  o m pone n t s  isolated. 

In  Chapter  5  'the  experimental  procedure  followed 
in  order  to  obtain  the  data  will  be  described. 

In  Chapter  6  the  different  methods  of  system 
identification  will  be  discussed  and  those  which  seem 
more  convenient  will  be  used  to  identify  the  motion 
parameters  of  the  mariner  class  hull  form. 

In  Chapter  ?  the  results  obtained  in  Chapter  6  will 
bo  discussed  and  in  Chapter  8  the  consequent  conclusions 
and  recommendations  will  be  pointed  out. 
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CHAPTER  2 

EQUATION  OF  COTTON 

In  order  to  develop  the  equations. of  motion  for  an 
ocean  vehicle,  it  is  convenient  first  to  define  two 
rectangular  right  handed  coordinate  systems ,  one  fixed 
in  the  surface  of  undisturbed  water  with  the  positive 
Xo  =  z   axis  pointing  downward  and  the  second  moving 
with  the  vehicle,  the  origin  of  such  coordinate  system 
being  arbitrary  (but  the  coordinates  (:',;',?„)  o^  the 
center  of  gravity  of  the  vehicle  in  this  system  being 
known)  pointing  the  positive  x-,  =  x  axis  to  the  bow  of 
the  vehicle  and  the  x0  -    z  pointing  downward  in  the 
equilibrium  position  (see  Figure  2,1) 


x?y 


x=z 


X,3X 


Figure  2.1 
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Th  e  orientation  of  the  moving  axis  system  with  respect 
to  the  fixed  axis  system  will  b"  expressed  "by  three  angles 
(the  "modified  Euler  angles")  defined  as  follows, 
assuming  that  initially  both  axis  systems  coincide  [3]s 

4"  -  angle  of  rotation  of  the  moving  axis  system 
about  the  z  =  z0  axis  (yaw  angle) 

9  -  rotation  of  the  moving  coordinate  system,  already 
yawed  about  its  y  axis  (pitch  angle) 

cp  -  rotation  of  the  moving  coordinate  system 
(already  yawed  and  pitched)  about  its  x  axis  (roll  angle). 

Defining  the  rotation  about  the  body  axis  as 
respectively  5^(about  the  body  x^  axis),  S<(sbout  the 
body  Xp  axis)  and  S^ (about  the  body  Xo  axis)  the  relation 
between  those  rotations  and  the  Euler  angles  will  be  {_3 J s 

$}     -  (j)no^cosf  -  9sin4- 

S     =  <£cos9sinvj~  +    ^cosf  (2.1) 

5^   =    -<j)sin9    +  *f" 

The  motion  of  a  rigii  body  in  six  degrees  of  freedom 
with  respect  to  on  axis  system  fixed  in  space,  when  the 
origin  is  a t  the  center  of  gravity  and  the  moving  axes 
are  the  principal  axes  of  inertia,  may  be  represented  by 
Newton's  equations,  namely: 


F  -  d    (muu) 
-   dt 


~       dt    "^ 


(2.2) 


where   m    is    a     iiagonal    matrix     ''.ose   elements   are    sc^ua]    to 


-15- 


th  e  mass  of  the  body,  u-  is  the  translational  velocity 

vector  of  the  center  of  gravity,  of  components  along  the 
moving  axes  u~  ,  vp  and  wp ,  F  is  the  total  vector  force 
acting  on  the  body  of  components  along  the  axes  X,  Y  and 
Z,  X  ^-s  a  diagonal  matrix  whose  elements  are  the  inertia, 
-Q?,    is  the  rotational  velocity  vector  about  those  axes 

Li 

of  components  p_ ,  q0  and  r^  and  M  is  the  total  moment 
acting  on  the  body  of  components  K,  M,  and  N.   Neglecting 
the  effects  due  to  rotation  of  the  ear+h  in  ocean 
vehicle  dynamics  Abkowitz  [_W~\   and  Tufts  [5]  develop 
equations  (2.2)  respectively  with  fixed  and  varying  mass 
o-n^  position  of  center  of  gravity  ^or  the  case  of  an 
arbitrary  origin,  the  first  using  vector  analysis  and  the 
second  Lagrange's  equations.   Therefore  only  the  resulting 
equations  will  be  presented  here  and  for  a  body  of  fixed 
"ass  and  inertia. 

X   =  mfu+qw-rv-x^Cq  +r   )+yr,(pq-r)  +  zr,(pr+q)J 

Li  Lr  [j 

2      2 
y  -   mTv+ru-pw-ypCr  +p    )+z    (qr-p)+x    (qp+r)j  (2,3) 

J  ■  j  Li 

Z   =   m[w+pv-qu-z_(p  +q    )+xr(rp-q)+y    (rq+p)  1 

""  G  Li  Li 

.    "        ■  2       ° 

K  =   I      p  +IYV(q-pr)    +    IY7(r+pq)    +    I      (q   -r    )    + 

(T      -T      )nr-inz    (v-DW+ru)    J.-   mv    (w-qu+pv) 

zz    yy;q      V  -g  ■  (2^) 

M  -    I      q   +1      (r-qp)    +    I^Cp+qr)    +    I      (r   -p    )    + 
(I      -1.      Vrn-mx^  (w-qu+pv)    +   mzr,(u-rv+qv/) 
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2   2 


N  =  if  +  I   (p-rq)  +  T   (q+rp)  +  I   (p  -q  )  + 

(T   -I   )pq-my  (u-rv+qw)  +  mx0(v-pw+ru) 

yy  xx      l-  ij 


(2.'0 


[u,A,u,^,X|Sf  input  | 


where  unsub-scripted  velocities  and  accelerations  refer 
to  an  arbitrary  origin  and  y.n,    y0  and  z_  are  the 

G     G  G 

coordinates  of  the  center  of  gravity. 

Note  that  F  and  M  will  in  general  be  functions  of 
position,  velocity  and  acceleration  of  the  vehicle  as  well 
of  properties  of  fluid,  geometry  of  the  body  and  of  the 
inputs.  For  a  body  of  given  geometry  and  mass  distribution 
and  in  a  given  fluid  then 
p   p 

v/1-  ere     u  =  [u  ,  v,w] 

-Q  -  [p,q,r] 

X  =  f~y ;   y  ,  Z  1 
—    -  o  '  "  O   O  • 

s  -  C  +  ,9.^] 

S^nco  the  purpose  of  this  investigation  was  the 
study  of  the  applicability  of  system  identification 
techniques  to  ocean  vehicle  dynamics  and  not  the 
identification  of  parar  ;ters  of  a  particular  vehicle 
(parameters  imbedded  in  F  and  M)  ;  t  was  decided  at  this 
stage  to  simplify  t^e  formulation  by  restraining  the 
model  to  one  or  at  most  two  5    grees  of  freedom.   It  will 


(2.5) 
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be  only  considered  heave  (translation  along  the  z  axis) 
and  roll  (rotation  about  the  x  axis)  with  zero  forward 
speed.   To  simplify  the  formulation  even  morr>  v„  ^nd  z 
were  ^ade  zero.   Therefore,  the  equations  of  motion 
for  this  restrained  mode  (yr,=  zr,=u=u=v=v=a>q=f=r=x  =v  =  ^^^) 

Lr    b  ~      x  0^0 

will  be 


m  w 


K  -  I   £ 

xx1 


and 


Z  =  Z[w-,w,z  ,p,pf  <(>  ,  input] 


(2.6) 


(2.7) 


K  =  K  [w ,  w ,  zQ ,  p ,  p ,  <j>  ,  in  rut  ] 
Looking  at  equations  (2.6  and  (2.7)  one  sees  that 
+  ^.~-r  ppg  almost  in  s  state  er^c^  form  sxce^t  "^v-^  f"i° 
fact  that  Z  and  K  are  ;,,'netionn  o^  w  and  o  and  that 


o 


=  d_  [s.J  =  d_  [4^o^9  cos^-9sinf*]   from  Eq .  (2.1) 


=  W  C  0  S  (\)  4 

o 


<p  +  v  3  i  n  cp 


see  Fig.  2.7. 


»i. 


Figure  2.2 
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(2.8) 


Since  H',  9  and  v  are  zero,  the  above  relations 
reduce  to 

<j>  -  P 

ZQ  =  w  COS  CD 

Ignoring  for  the  moment  the  presence  of  w  and  p  as 

arguments  of  Z  and  K  the  cor-plete  set  of  equations  in 

state  space  form  would  be: 

w  =  1   Z(wfw,zo,p,p,(j>f  input) 
m 


4> 


z   =  v.*  cor 
o 

p  =  1   K(w,w,zn,p,p,£  input) 

XX 


(2.9) 


o 


4> 


=  D 


In  order- to  be  able  to  proceed  and  apply  the  existing 
identification  techniques  to  ocean  vehicles  (specifically 
surface  ships),  it  is  necessary  to  isolate  w  and  -that 
the  dependence   of  Z  and  K  on  w  and  c  is  linear.   Therefore, 

A  J-  * 

the  next  step  will  be  to  assume  the  existence  of  a 

Taylor's  expansion  about  an  equilibrium  po.int,  to  be 
taken  as  w=w= zq  —   0   of  the  form: 


Z  -  z(o)  +  SjZ 
6w 


6Z. 

5  inn 

3 


6Z 

"o  -1-  6^ 


n  + + 


(2.10) 


w2+ 
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I  ill 
2!  6W6w 


2 

ww  +  I  6,'z 
2!  6w6z 


1   6?Z 
wzo  +  2!  6*6* 


V.'D+ + 


(2.10) 


and  in  the  same  way  for  K, 

It  is  known  from  potential  theory  [5 1  that  for  a 
body  moving  in  an  infinite  fluid  away  from  the  surface, 
there  are  no  dependence  of  force  or  momenta  on  products 
of  accelerations  with  velocities   or  accelerations. 
It  will  be  necessary  to  accept  this  assumption  for 
surface  vehicles  without  any  other  justification,  and 
neglect  also  all  other  terms  where  acceleration  is 
multiplying  any  other  variable. 

With  the  vehicle  in  equilibrium  also  7(o)  =  0, 
Further  simplification  can.  be  achieved  if  we 
cans  id er  each  element  of  the  expansion,  which  will  be 
considered  in  the  next  chapter.   It  will  be  enough  to 
say  for  the  moment  that  due  to  port  and  starboard  symmetry 
(both  geometrically  and  inertial)   the  vertical  motion 
of  the  vehicle  does  not  produce  any  rolling  moment,  at 
least  to  first  order,  and  therefore  those  terms  car  be 
drcoped  from  the  equations  further  simplifying  them  '  '-•■  ' , 
Similarly,  and  to  first  order  only,  rolling  motions  will 
not  produce  any  vertical  force  or,  in  other  words: 


6Z 
6- 


6Z 


6f 


6K 
6w 


6K 

6  w 


6K 
67, 


0 


(2.11) 
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Also  the  square  terms  62Z  and  65K  will  be  nonexistent. 


T" 


6W6      6p* 

Instead,  consideration  must  n°  given  to  frictional  and 
vortical  effects  which  are  of  the  form  c-^wlwland  c  p  ]p|  . 
Although  these  terms  are  not  analytic,  they  will  "be 
included  so  that  the  model  represents  a  closer  approxi- 
mation to  the  vehicle  physical  behaviour. 

The  form  of  the  equations  will  then  be: 


w  -«*- •  u 


w+Z 


w|w| 


+ input 
w  cos  CD 


r  t^     n    m 

w  +2,   z  +...  +  _  1    Z„  xn  (x.)  (x.)  -*■ 

z    o 
o 


(rw-n)!   Xi  J    x 


(2.12) 


1 


-K'+I 

P   xx. 

+ input 


K    T3+K     n  i 
'V   Plpl 


-  _1    K  n  m  (x.)n(x.) ' 

( m+n )  ;  y-    X.    J 

J  1 


d>  = 


=  P 


where 


zx  =  & 

3 


and  K 


x.         6x 
In  the  next  chapter  each  term  in  the  a"bove  equations 


which  is  of  first  order,  quadratic  or  at  most  cubic  will 
be  discussed  in  relation  tc  its  nature,  order  of  magnitude 

and  methods  of  predicting  its  value. 

In  summary,  from  Newton's  equations  and  for  a  case 
of  "restricted  motions,  assuming  the  existence  of  a  Taylor 
series  for  the  force  and  moment  and  linear  dependence  of 
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foroe  and  moment  on  acceleration  it  was  possible  to  obtain 
a  set  of  equations  in  a  state  space  representation  to 
which  some  system  identification  techniques  may  "be 
applicable. 
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CHAPTER  3 

r-'OTTO?1!  PARAMETERS 

For  convenience  of  presentationi  the  ^c.pi;s?ior  on 
+h?  motion  parameters  will  be  subdivided  into  three  main 
groups,  namely: 

i)  motion  parameters  resultant  of  differentiation 
of  force  or  moment  on  the  body  with  resoect  to  position 
(7  ,d>)  only  (hydrostatic  coefficients) 

ii)  linear  terms  resultant  from  differentiation  of 
force  or  moment  on  the  body  with  respect  to  acceleration 
and  velocity  (added  mass  or  added  moment  of  inertia  and 
d;  mping  coefficients) 

iii  )    non~l  inear  hvdrodvnannc  fe~^~~F 
3«1~  Hydrostatic  Coefficients 

These  coefficients  are  dependent  only  on  the  geometry 
and  mass  distribution  of  the  body,  for  example,  shape 
and  relative  oosition  of  the  metacenter,  center  of 


■k 


DU 


oyancy  and  center  of  gravity,  for  the  restoring  noments . 


The  coefficients  •  ^"",  calculated  1  ;  3  on  the  lines 
plan  of  the  hull  for-1. , 

With  the  ship  in  the  upri  ht  position,  the  buoy;   t 

force  B  (z  ,0)  is  calculated  for  each  waterline  (each 

o 
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selected  value  of  z  ).   Then  at  each  water-line  an  angle 

of  roll  is  simulated  in  the  body  plan,  being  the  point 

of  rotation  the  assumed  position  of  the  center  of  gravity, 

and  the  resultant  buoyant  force   3  (  z  ,  <f>  )  and  the  righting 

arm  GZ(zQ,<p)  are  calculated  (for  details  see,  for  example, 

"Principles  of  Naval  Architecture"  edited  by  J.  P.  Corastock, 

Chapter  2).   The  righting  moment  (  -K(z  ,<p)  )  is  then 

given  by  the  product  of  B(z  s<f>)  •  GZ(z0,^>).   The  vertical  force, 

positive  downwards,  is  the  difference  between  the  weight 

of  the  body  and  the  buoyant  force,  or  in  other  words, 

considering  the  component  of  the  force  along  the  body 

v»  ay  is  ,   .' s  haves 

Z-(zo,<|>)  =  (W  -  B(z0,£))  cos  <j> 

K(z0,f)  =  -B(z0,<f>)  •  GZ(z0,cj>) 
where  W  is  the  weight  of  the  vehicle. 

The  values  of  Z(zn,a>)  and  K(znt<f>)    are  usually 
presented  tabulated  or  in  graphical  form  (see  Appendix  1), 
In  order  to  get  the  hydrostatic  coefficients,  lerivatives 
of  (3.1),  it  is  necessary  to  obtain  a  mathematical 
representation  of  those  surfaces,  for  example,  a  polynomial 
fitting  by  the  method  of  least  squares.   One  such  scheme 
was  developed  by  Kerwir  £6~]  and  I  rplemented  by  Parlssis  £?~l 
using  Legendre  Polynomials  to  represent  ship  forms  which 
is  equa 1 3  y  app] icable  to  ( 3  .3 ) . 
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If   or5    is    interested    in   using    in   the   state    equations 
just    the   linear   terms    such   a   procedure    is   not   neer'--r-nry. 
It    is    easily   seen    that 


6Z(z0,j>) 
6* 


,4-o 


=       6K(zol<f) 


5z 


=    0 


(3.2) 


Zo=<M 


due  to  port-starboard  symmetry  of  the  hull  form  and  that 
6Z(z0,<f>) 

o=<M 


6z 


TgA'VL 


6K(z  ,<{>) 


(3.3) 


-W  •  GM 


cj>-o 


where  A-,rT  represents  the  area  of  the  waterolane  at  the 
design  v::"  t^rl  ire  (zo=i0)  are  GM  is  the   9t<  centric  height 
(distance  from  the  metacenter  to  the  center  of  gravity). 

Due  also  to  port  and  starboard  symmetry. in  the  non- 


linear formulation  the  term  1    6  K(z  ,  ) 

21  A   2 


IS 


zo=<p=0 


equal  to  zero.   Therefore  in  "''his  case  the  cubic  term 
will  be  maintained  in  the  equations. 

3 .2    -    ArM  ~d  Mass  ,    A^d-^  foment  of  Inertia  dnd  Da-nirm, 
Co-" f  f  5  c i  e'nts 


From  the  eig^t  fT^-s  that  form  this  subgroup  four 

wil]  be  zero  again  due  to  symmetry  of  the  vehicle.   Those 

wi.ll  be  the  cross  coupling  +  ?rms 

-Z*=Z  =-K«=-Kw  -0  (3-lL) 

p   p   w   w 

and  the  terms  to  be  retained: 
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•Zw 

(3.5) 

•KP 


So  it  is  seen  that  when  only  linear  terms  are 
considered  the  equations  of  heave  and  roll  are  completely 
decoupled  and  may  be  studied  independently.   If  the  six 
degrees  of  freedom  were  considered,  the  existence  of 
port  and  starboard  symmetry  would  cause,  when  retaining 
only  the  linear  terns  in  the  equations  of  motion,  the 
decoupling  of  the  equations  into  two  groups  ["3<'n* 
coupled  surge  (translation  along  x^  axis),  heave 
(translation  along  the  x.o)  and  pitch  (rotation  about  the 
Xp  axis)  -  the  so-called  motions  in  the  vertical  plane, 
and  coupled  sway  (translation  along  the  X?  axis),  roll 
(rotation  about  the  Xi  axis)  and  yaw  (rotation  a/bout 
4"^e  Xo  axis)  in  the  horizontal  plane. 

Added  mass  and  damping  coefficients  are  in  general 
not  known  except  for  some,  special  types  of  idealized 
forms  and  when  the  vehic] es  are  excited  into  some  special 
types  of  motions.   One  such  type  that  have  been  studied 
quite  extensively,  especially  in  the  domain  of  surface 
vehicles  or  vehicles  near  the  surface,  is  that  of  a 
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vehicle  excited  "by  waves  into  an  oscillatory  motion.   t^10 
path  usually  followed  [3.8,9]  is  to  assume  the  fluid  to 
be  inviscid,  incompressible  and  irrotational ,  therefore 
assuring  the  existence  of  a  velocity  potential  that 
satisfies  both  the  Laplace  equation  and  the  Bernoulli 
equation.   Considering  the  motion  harmonic,  the  solution 
is  sought  as  the  product  Re  :^f (  Xn  ,x9  ,Xo  )  e     '  where  <*>   is 
the  frequency  of  oscillation.   Considering  the  linearized 
problem,  t(x-i  tX-2* x3 )  ^s  assumed  to  be  the  sum  of  eight 
other  potentials  being  the  total  motion  the  superposition 
of  each  due  to  each  of  the  potentials:   two  of  them  due 
k0  khe  ■  ction  of  the  incid  ~^J~    wave  plus  !  iffr;  :  t  ;  :  vave 
on  a  fixed  vehicle  and  six  due  to  forced  oscillations 
of  the  vehicle  on  calm  water.   The  fluid  is  assumed  in 
all  cases  to  extend  to  ±°°  in  the  X]_  and  x?  directions  and 
between  Xo=0  and  +°°  or  between  x^~o  and  x^1,  for  the 
finite  depth  case.   The  linearized  bo u  idary  conditions 
will  be: 

a)  On  the  surface  Xo=o,  the  velocity  of  the  fluid 

•  5  -  _£  where  $   is  th  =  surface  slevation,f or  all 
6X3    6t 

potentials  and  is  ;qual  to  -^    ~P.    u:  ing  the  line;  -;"  d 

a-  '  ■ 

o 

Be rn  ou .1 1  i  eq  nation. 

b)  On  the  body  for  ;  h  ;  incid  ;n1  and  diffracted 
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potentials   j^£  =   Q  =    6 (*rw+ yd )    v/here   n    is    the   unit   normal, 

6n  5n 

For   the   six    other  potentials    the   norma]    component 

of   the   velocity  will    eQ.ua]    the    velocity   of   the   body    in 


that    direction    or 

R 


=R  -  '^ 


.ia 


e 


-i^t 


f.i 


i"t 


6n 


-    f -s. 

■1    3 


or  being 


(s.      is    the  arnolitude    of   the    oscillati 


on 


in    the    j    direction)    h^- 


6n 


=    -i<*>f  .s  • 

t  "  na 


(3.6) 


•/.'here   f.    (  j   —    1 ,  ....  6 )    are   the   generalized    direction   cosines 
of   the   body   surface 

fi  (j=l,2,3)  -  cos  (n,x.) 


f 


U 


x2-3  T2 

=  x  f  -Xif 

3  1  ]  3 


f6 

X    f     -X    f 

12      2   1 

c)    At 

X~     =     +e» 

V$=   0 

x,,   -  +h 

t?  "   o 

(3.7) 


2    2  i 


d)  At  r  =  (x-i  +x   ) ":    the  radiation  condition  must 
be  satisfied  for  the  six  potentials  and  for  lif fraction 

potential  nv : 


lim  r" 
r-**> 


^    ■     2  0) 

g 


6r 


0 


j  =  1...6,d   (3.8) 
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M  ;o  Prom  Bernoulli  equation  the  hydro'dynamic  pressure 

w  ill    "be 


"0    = 


T3^ 


-  i  w>-i- 


(3.9) 


The  hvdrodvnamic    force    in    the   k   direction  will   b< 

6 


Fk   =   -H     pfkd.s=Re 
o 


-  iiopt 


<Jt 


s^r  «-^v 


(3.10) 


The    total   hydrodynamic    fore?    o^    the   body    in  a   given 
direction  k  was   defined    to   he    in  Chanter   2   for  .^otions 


n    r^ode    i:    F,    .    =    6Fi,  •      s.    +    6Fv;      s 


k.i 


v  i 


J 


ki 


J 


6sj 


6s 


;  P  1  r\  ,v 


J 


Ref-i<os  • 


■i<ot 


a 


] 


and   changing  the  nomenclature   to 


«^_ 


6s  . 


-Akj. 


-    -p. 


6s 


"k.i 


-  5  o>  t. 


3 


Ro[~(^An   .   +   i  «   B,  -)s.      e  (3.3  "0 


and  the  total  force  in  the  k  direction  was 

6 
P,  =  .s,  F,  •  +  P, 

Fv. ,  being  the  excitation  in  the  k  direction  (the  hydro- 
static  "^orce  was  excluded  from  the  fores^oii  oression 
as  it  was  from  the  Bernoulli  equation  (3.°)  ). 
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Comparing  the  foregoing  expression  with  (3.10),  the 


following  is  obtained: 

_  "  03*- 


wk 


o 


(o  A,  .  +   ioB,  . )s  .   =  - 


K.I 


v-ja 


Is  TJ 


f,ds 


(3.12) 


(3.13) 


and    therefore  we   obtain 


Ki 


s  .»  u>     S 

ia  o 


s  .v       S 

ja         o 


'      LIjJ  k 


(3.1*0 


that  would  give  the  value  of  the  hydrodynamic  coefficients 
if  the  velocity  potential  was  known.   Since  that  is  not 
the  case,  it  is  not  possible  to  use  the  above  -- ;  aticiis. 
Instead  one  looks  for  possible  relationships  between  the 
hydrodynamic  coefficients  to  get  a  better  insight  on  the 
s olu t i on  o f  the  pr o bl em , 

For  example  using  the  second  form  of  Green's  theorem 
to  a  velocity  potential  «P.  end  to  another  velocity 
potential  ^P  with  k  ^  j  we  obtain  since  botl   x  '  itials 
satisfy  Laplace's  equation:  [3] 


A 


k 


A 


J* 


% 


B 


(3.15) 


Also  from  [°"1  and  [lOj  the  following  expression  is 
obtained  relating  =  Ided  mass  and  danroing  coefficients: 
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OO 


A1  .(u>)  -  Ak.(-)  =2|  [Byi(  f  )-B,.(oo)]   d£ 


kj 


OO 


(3.16) 


(-)  =  1  X  U.  .(f  )-A..(oo)]1fdI 


b  .  ( co )  -  Bv  • 


Maybe  the  most  Important  relations  are  the  so-called 
Haskind  relations  |_ll[]  based  also  on  the  second  form  of 
Green's  theorem  when  applied  to  equation  (3.12)  after  the 


ixpression  for  f   in  (J, 6)    or: 


f.   -   -A^V 


-6  Tk 


5n 


ka 


is    substituted    in    (3.12): 

_1-t  XX      (fw  +fa)6^k  as 


T? 


=    k  e 


'  Y>  i~>  p  v, 


wk 


f-  e 


s 
L  ka 


6n 


(3.6) 


(3.12) 


fO] 


/J*/    (f&M-<R*f)    dr   -  <3        (J>£f;      -     4.  tV?) 

R  s  '6n  6n 

o 

if.  fd 


(3.13) 


Since  "both  potentials  satisfy  Laplace's  equation  the 
left  hand  side  is  zero  and  therefore: 


SI 


vO  _   &fk 


6n 


U  il 


£fe  ds 

6n 


(3.1*0 


nd  from  the  boundary  condition  on  the  body 
6fi    =   -  6*fw 


(3.15) 


6n  6n 

The  following  wil]  result  for  the  excitation  force: 
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■pi         _  n  - 

rwk 


p     e"iu>t    r;   [<fy   5ft   -  fk    sfwl   ds 


(3.16) 


eliminating  the  diffraction  potential. 

Following  [ill  ^or  2-dimensional  bodies  (under  "beam 
seas)  the  following  relations  are  obtained,  "based  on 
(3.16): 


p, 


kk 


r 


Pi 


and  therefore 


(3.1?) 


Pv-  =  *  [_fdBkVl 


to 


kk 


k=   2,3,4 


where  />a  is  the  amplitude  of  the  incoming  waves  and  Fv 
is  the  amplitude  of  the  exciting  force.   Being  the 
darrein*?  coefficient  calculated  fron  the  energv  radiated 
-  +  infinity  it's  possible  to  calculate  the  magnitude 
of  the  exciting  force.   However,  it's  phase  is  not 
determined . 

From  the  foregoing  it  is  seen  that  the  analytic 
solution  for  the  general  case  of  a  three  dimensional  body 
does  not  exist.   Only  the  cases  of  the  sp:  bt   s  end 
Qllxpsoi  Is  have  been  solved  [l2J.   It  would  be  possible 
to  solve  a  problem  numerically  by  such  a   rocedure  as  the 
finite  element  technique  [3!,  hut  nothing  is  available  yet. 

In  order  to  get  information  for  practical  applications- 
it  is  necessary  then  to  work  into  two  dii   isional  problems, 
The  strip  theory  for  ship  motions,  introduced  by  Korvin- 
Kroukovsky  [l3i  1''  »  15l  makes  use  of  two  dimensional 
solutions  for  each  transverse  section  o^  the  body 
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detornining  the  characteristics  or.  the  flow  and  conse- 
quently the  local  hydrodynamic  coefficients  and  then  the 
results  are  integrated  over  the  length  of  the  "body  to 
obtain  the  hydrodynamic  coefficients  for  the  whole 
vehicle,  applying  then  a  correction  factor  to  take  into 
account  the  existence  of  a  three  dimensional  flow. 

The  solution  in  two  dimensions  has  been  obtained 
for  the  three  possible  modes  of  motion:  sway,  heave  and 
roll.  One   way  to  obtain  the  solution  is,  using  conformal 
mapping  to  transform  the  cross  section  of  arbitrary  shape 
of  a  two  dimensional  cylinder  into  the  unit  circle, 
j-q   assume  the  velocity  potential  to  be  s  sum  ^"f  mu]  bipole 
potentials  whose  coefficients  arc  a ,5  justed  to  satisfy 
the  boundary  conditions  on  the   cylinder  [_l6j.   For 
example  for  heave  the  multipole  expansion  will  be  formed 
by  the  sun  of  a  source  potential  at  the  origin  plus  the 
multipole  potentials  while  for  sway  and  roll  instead  of 
a  source  there  will  he  a  dipole  since  the  heave  potenti;  ] 
is  svmmetrical  end  the  other  two  asymmetrical « 

,T1he  conformal  transformation  will  be  of  the  Form 

rsi 


a 


-  f 


N 

+    2 

u-0 


a, 


'2n+l 


~(2n+] ) 
f  (3.18) 
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which  maps  an  arbi trary  section  in  the  z  plane  into  the 
unit  circle  .in  the  ^plane.   From  the  section  offsets  it 
is  necessary  to  find  the  coefficients  &_,•,,   A  method 
of  determining  those  coefficients  was  developed  "by 
Bermejo  [l?]  using  Theodorsen's  method. 

Another  way  of  solving  the  two  dimensional  problem 
is  to  represent  the  velocity  potential  by  a  distribution 
of  sources  on  the  contour  of  the  section  [l?1.   The 
intensity  of  the  sources  is  determined  applying  the 
kinematic  boundary  condition  to  the  underwater  portion 
of  the  section.   Once  the  potential  is  obtained  the 
oressure  is  calculated  usin0-  equation  (3»9)  —  linearized 
Bernoulli  equation  -  and  the  force  is  obtained  by 
integration  from  where  the  added  mass  and  damping 
coefficient  are ' calculated  (En.  3.9  to  3*12), 

Based  on   this  method  a  computer  program  was 
developed  [19]  which  performs  the  above  calculations 

5  was  used  in  this  investigation  for  the  computation 
of  the  initial  conditions  of  those  motion  parameters,, 
Listings  of  those  programs  are  presented  in  A]    iix  ?. 

Still  another  way  of  calculating  hydrodynamic 
coefficients  is  to  use  approximate  formulas,  based 
some  of  them  on  the  same  principles  as  exposed  above, 


-34- 


and  some  in  model  experiments  and  full  scale  information. 
Some  of  those  formulas  are  presented  by  Blagoveshchensky 
[20]  and  v/ere  used  here  to  calculate  A^,  A^;,  ,  B^o  and 
^UU   anr^  "^e  resul"ts  compared  with  these  obtained  with 
the  computer  program  mentioned,  believed  to  be  more  exact. 
Those  results  were  also  compared  with  the  ones  obtained  by 
Hines  and   Gies  [21j  in  experiments  with  another  model 
of  the  same  hull  form  ,   Those  results  are  presented  in 
Appendix  3. 

3.3  ~  Non-linear  Hydrodynamic  Coefficients 

In  this  category  are  the  terms  of  second  order  in 
velocity  and  cross  terms  in  velocity  and  position. 

Information  en  these  types  of  parameters  is 
practically  nonexistent.   Only  two  references  to 
factional  damc^n?  in  rol]  were  found  [ 20 ,  22*1.   Since 
that  is  the  ease  where  frictional  and  separation 
resistance  has  a  bigger  effect  when  compared  with  the 
other  modes  of  motion,  the  expressions  suggested  by 
Blagoveshchensky  was  used  since  it  was  directly  applicable 
to  ocean  vehicles. 

The  expression  given  by  reference  '  22  '  relates  th^ 
friction  coefficient  of  oscillating  cyclinders  to  an 
equivalent  Reynolds  number: 


cf  =  1.328 


T  OO     w*  A 


m 


<"P 


.5 


(3.19) 
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where      v-     -   radius    of   the      cylinder 

9       -   rne^n   amplitude    of   the    oscillation 

m 

T   -  period  of  oscillation 
v  -  kinematic  viscosity 
This  formula  is  the  same  as  blasius  solution  for 
C+  in  laminar  flow  if  the  Reynolds  number  is  taken  ?,s 
Rey  =  3.22  r-5  8m  (3.20) 

To  obtain  the  drag  coefficient  for  heaving 
oscillations  a  similar  method  may  be  used,  although  with 
rougher  approximations  since  the  vehicle  is  not  a  two 
dimensional  body  and  the  sections  do  not  resemble  circles. 

He^e  a  "strip  theory"  approach  v.";,'7  used.   Results 
were  sought  for  each  station  and  integrated  over  the 
length  of  the  body. 

For  each  station  it  is  assumed  that  the  resistance 
is  due,  '/.'hen  moving  into  the  water,  entirely  to  skin 
friction  end  when  moving  away  from  the  water,  entirely 
to  separation.   Therefore  for  the  first  case  the 
diameter  of  the  equivalent  cylinder  will  be  d  =  i/tt 
where  1  is  the  length  of  the  underwater  part  of  the 
contour  of  the  section,   For  the  second  case  the  l.ia 
of  the  equivalent  cylinder  will  be  equa]  to  the  1 
of  the  section.   The  velocity  used  to  calculat  ■  the 
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Reynolds  number  will  be  the  "root  moan  square"  value 

based  on  the  tests  perf ormc^ .   Vol!]'"'"  of  skin-friction 
coefficient  and  form-drag  coefficients  are  obtained 
from  reference  [23^.   Results  of  this  procedure  together 
with  the  results  of  the  coefficients  as  calculated  by 
the  method  indicated  in  reference  [  2QJ  are  presented  in 
Appendix  3« 

With  respect  to  the  other  coefficients,  the 
assumption  may  be  made  that  they  arc  indeed  very  smell 
and  consequently  the  assigned  value  to  be  used  in  the 
identification  procedure  is  zero. 

Summarizing  in  this  chanter  the  si  i    Dt.io"    ra— 
- ,-.-'-  gv>g  w°r  2  c  ''~''  s  ]  d  e^od   t]  g  i  r  ns  bu^e   "  '  ■    ■   s  i  d 
the  available  methods  of  prediction  outlined,   '^v  e  results 
obtained  for  the  chosen  hull  form  are  presei  bed  in 
Appendix  1,  2  and  3« 

In  the  next  chapter  the  exciting  force  and  moment 
will  be  considered  thus  completing  the  description  of 
the  mathematical  model. 
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CHAPTSR  ll 


E/CTTATTCN  FRO!'-  T::E  E"VIRO?':r.,Tgr-:T 

Fro~  all  possible  tyros  of  excitation,  only  forces 
and  rodents  due  to  the  action  of  two  ■'     "   a]  plane 
progressive  waves  "re  considered  because  (1)  t!  =y  ere 
well  defined  in  mathematical  terms  and  (?)  easy  to 
generate  in  a  towing  tank  experiment. 

^h r    expression  for- the  velocitv  potential  is 
obtained  "Pollowiner  a  oath  similar  to  the  one  followed  in 
C,iort?r  3  for  the  calculation  of  relatj  c  is  v  tv/een  the 
linear  hydrodynamic  coefficient?;',  namely  that  the 
veloc  ity  pot ent ia]  sat isf  i    s  botr.  Laolacci      Bern ou 11 i 
equations.   Taking  the  case  of  ar  ;   'inite  fluid  in  all 
directions  thn  folic/;'"  "  jxores  >ions  hold  r 1  c '.  2?i! 


s/<£  w 


-  0 


d(5-x  ) 


6t      6x 


*3=f 


>~* 


=0=-  ^  4>-  i  '  (v<|w)3  -*•*> 


{,4wl) 


3 


x  =0 

3 
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<£v:  =  fin  it; 


■y  — ^  oo 
y'3 


^rifTT,' 


V+.D 


where  ^  is  the  surface  displacement. 

The  solution  is  sought  by  separation  of  variables 

and  assuming  a   harmonic  dependence  on  time,  or  in  other 

words,  and  for   two  dimensional  behaviour: 

-itft-i 

5 


W   ( ; 


,t)  =  Re[<f2(xs)  .pr.3) 


v    v-  „   y  . 


C^.2) 


Plugging  this  solution  back  into  Laplace's  equation 

we  obtain 


=  -k' 


C^.3) 


D  0  S  ] 


To  satisfy  the  last  boundary  condition  k3  must  be 
tive  being  the  solution  for  [  (v2)  anc^   j  (x  ): 

ikx. 
f?tx2)   -  C   e    ' 


-1  KX. 


+  C   e 

2 


kx  -kx 

:   3    +  D   e    3 

2 


lf(x  )   -  D   € 

[3    3  1 

1  n d  to  s a  t  is f y   v1  CD  =  0    x 


C*.*0 


3  1 


=  0 


f  (x  ) 

|3    3 


=  D 


■*-3 


and  therefore 


t,(x 


■Vy. 


,x  ,t)=Re[  e   '  3(C  e 

u    •    1 


(kx2-«t)    -i(kx2+«*t) 


+C 


2    3 


) 1   (V5) 
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Th   olution  with  C  _  0  corresponds  to  waves 

progressing  in  the  oositiv0  x   direction  and  with  C   =  0 

2  1 

to  waves  progressing  in  the  negative  Xp  direction. 

Assuming  the  later  case  and  the  existence  of  an 
arbitrary  phase  angle  §  to  be  determined  later  from 


the  initial  conditions: 

I  ,  +>  _   „   -[kx  +  i(kx  +«Jt  +6  )] 

q>  (y  ,x  ,t)=Re   C  e  u   3        2  J 

J  w   2   3 


(*K6) 


Combining  the  two  boundary  conditions  on  the  surface 


(x  a   0) 
3 

e   6t^~ 


1 


CO 


5  w 
6x 


k 


x  =  0 
3 


[^  .?) 


which  determines   t!       value   of  k,    the  wave  ?r. 

Us  inf 


J  6t 


0 


i  fkXa+^'t    +  6   ) 
=    1      Re    f    ~'.W  C    e 


«Kl 


w. 


'       (kx  -!-   uit+6) 

2 


and    for   the    form 

-^  =  -&      sin(kx  ■f«t+ 6    ) 
J  3 


(fc."8) 


,..'_  ?re   y      •..    ^g   arpp]  itude    of   the   surface   disp]      :e         :  ■ 

C      -      rj_  ('''   .9) 


OJ 
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an< 


$  (x    ,x    ,t)=R< 


n 


(/.  (i  0) 


from   the   Haskind    relations,    Eq .     (3.l6) 


F     =  R  < 


6r. 


where 


L   Cp~  o 

vf(x1  ,x    ,x    )=     $w(x,  ,x    ,x    ,t) 

Iw  2  3  ■*•••-  2  3 


on 


i^t 


fpvifo    first   term   in   the    in  t :  irrs  ]  "^~      Fc  " 

exerted    by  the   pressure    in   the  wave   withoul     taking  ; 

of   the   disturbance   "by  the  presence   of  the   ship   and    is    !'•  e 


only   excitation  accounted    ° o r 


>v 


•-  -    Fr  ou  d  s  -  v  r i  i  o  f' : 


hypothesis.      The   second    tern    is    the   contribute  )n    of   tl  = 


r-  "*-  1    i"V 


rp0    coFiuute   the    first   term   in   the    exci  '  '         f  rce 
(the    Froude-Krilof f   force)    first   replace 

6  fi     ^  i :Uf  .  s  .  (3.6) 


:Vf  -s  . 


(F    .)..     =R. 
w.i    fk 


u)  f 

i 

Pg^-SS       -~      3sin(kx   +t»>t+6)f.ds 

'  Q  2  u 


s 

o 


('•-ID 


The    second    term    (diffractior     force)    ii      •  ■-.■• 

noting    thai: 

,o  „        -  \  kx   +i  ( kx    +6)1, 

5Ty;  =     -w-fa         -     3  a      '    (f  +if    ) 


6n 


v:  t  '  d 


F7~ 

ia 

V 


-  ]    ^  t  n  r> 


?, 


1     ■.-       4.'    (W 
•-""o-  "2 


(f  +if   )ds 

2  3 


» 

(a. 13) 
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°nd 


F  .  =  (F   ) 

w  i   v  w  i  f  k 


(Fwj)d  (".!') 

Therefore  "to  calculate  "the  bote  ]  excitation  force 
the  results  from  the  sub-problem  of  forced  oscillatio  ■ 
of  the  body  in  calm  water  must  be  known,  namely 
T-(x-,  ,x  ,x  ).   What  was  said  about  the  way  of  calculating 

U    x    2    3 

oQ^o.-j  mass  end  damcine,  coefficients  3 r'   acol  Lcable  hers. 
Besides  it  is  "possible  to  calculate  this  second  term  In 
tve  exciting  force  at  the  same  time  as  the  added  mass 
and  damping  coefficients  are  calculated,  since  "^he 
potential  is  also  obtained  in  the  rode  of  interest.  Or, 
5r»(js  }  v.'o  ol   in  ^"o"^  each  its  b?  e%~ 


1^  Q- 


P  F    .  K 


Qo 


OA    i   >y 


P 


,,      -iwt    p  -[kx3+i(kx3+6)l 

u)  .  e 


(f  +if   >d] 

2  3 


[Il.l5) 


where  C      is    the   contour   o_r  the   station    and    d]     is  'the 
o 

I  9]np"th    element   a]  one   that   contour. 

These   two    :"    ensiona]    res  ilts    for    each    s tat'?. on   "re 
then    integrs  !"  ~  i    over    j~''.e   length    o  IT    bhe    vehicle. 

Al thou°"h     th is  1     '  >    to  be    tl  i  r      >ui1 

for   the    computation    of         '      xcit?         fore       ,    that  :" 

used    in   this    invest  i^  ition  due   to    bl       f  !  at    the 

c omputation    of    ?ouat  ion    ( ■''  .  3  5 N     v " r^    noi     included    in    bhe 
computer   scheme  used    in    tl       :a]       I     bion    of    ;ectii     "! 
added    mass    and     1;     pi  Fficients.       Cns1        ;,     bhe 
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procedure  presented  by  B] '  c)   ei  sky  [?nl  was  used 

since,  although  less  correct,  r<   r  lered  the  computations 
much  easier  to  perform. 

First  it  is  assumed  that  $£_   «1,   B/^«l,  ?/;\  «1 


where  ^fil  is  the  wave  amplitude,  ^  ; s  the  wave  length 
A~  2-rr/k,  B  and  T  are  the  beam  and  the  draft  of  the 
vehicle. 

Next,  it  is  assumed  that  the  principal  part  of  the 
excitation  is  due  to  change  in  buovancv  due  4~o  the 
presence  of  the  wave.   Then  some  multiplicative  correction 
c ^ ?f f ic i en ts  are  calculated  to  account  for  the  curvature 
of  wave  (fin: be  beam  correction)  and  ror   th       1  nee 
of  the  h^drodynamic  oressure  on  x   ( Finite  draft 

3 

correction).      Each   correction    is    calculated    separately 
rirfi    assumed    independent    from   each    other    (for   example, 
when   calculating  the   beam   correction,    T/\«]    and   vice 
versa) . 

Under  these   assumptions    the    linear   equations    ~"^~r 
h  e a v e   and    roll   w ill   t a  k e   1; h  e   f or m : 

a)    heave 
w   =  Z  w         f  \z  (■:       -  fc) 

..'7   •  J-m 

(k.l6) 


■Z'+m 
to 


■Z'+m 
w 


w 
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wher^  jc  is  "the  surface  displacement  corrected  for  finite 
beam  and  draft  effects 

^c   «*    zb      zt 
=   a  sin  (ky  +   t  +6) 

b)  roll 


o  = 


"D 


-K'+I 

O   XX 


•[  K  ^  -  W.BG](4~  5l  '  +   W,EC   4 


(^-17) 


d)  -  p 


iirVi  i  V  f 


-K'  +  I 
p   xx 


-K'+I 

P   xx 


v^  io  !_,>"'"  wave  slooe  corrected  for  fir  5  ^~    bean  and 

d  ra  f t  g  ^"  "^  °  c  t  s 

V  =  5;;      C0b  C^t 

JP    -  if   -  £  CO  ;(Vy  +  »t   +6   ) 


C   ,  ,    C      '    C      ,    C,,    are    the   nultiplicati"1  :-   corrections 
zb        zt       p>v      ^t 

given   by  [20  ' 


Zt} 


1-   1.73  C      (e/a)2 


G        -i-o    (kT) 
zt 


q  f1^)3.        c  (kT)3 

2(2-C    )  ^T^H")" 


(a. 18) 


Vt 


6    G    3 

(i+c  )(i+2Cmr 


-1 


(kT    )+         i  .  5   G    3 

m 

(2~G   )(2+e   ) 


(kT    ) 


Cjt  (kT)3 

3T>    :m)(3-c]T1) 


wb  ?re 


C  is  the  vvaterlin.e    coeffici  -  *    G      =  A       /-T 
w  w 

G  is  J  '•        ndshii]     i  set  ion    coefficient   G      =    ^       /BT 

k  is  the  •    :  -' "    mi    ber .    • 
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Th  e  details  of  the  calculation  of  these  correction 

coefficients  may  he  found  in  Reference  L?0J  paragraphs 

27  to  31. 

In  Equation  *K17  the  term  ','/»  PI       represents 

p   xx 
the  contribution  due  to  the  inclination  of  the  ship  when 

Q-   y     since  the  weight  of  the  vehicle  and  its  buoyancy 

will  not  be  in  the  same  vertical  thus  causing  a  heeling 

moment  W.PC-   sir.  d>    being  W.BG-<£    the   first    tpr^.    in   the 

Taylor's    expansion.        W    is    the   weight    of   the   vehicle   and 

gQ    is    t1"! °   distance   from   the   center   of    '"< — ivitv   J~o   the 

center  of  buoyancy. 

Although  Vc. e  expressions  above  for  the  exciting  force 

and  moment  are  not  as  correct  as  Equations  k-.ll,  ^.13 

and  ij-.l^,  the  simplicity  of  their  calculation  decided 

their  introduction  in  the  equations  of  motion, 

When  considering  the  non-linearities ,  if  the 

exciting  force  can  still  be  approximated  r,s  expressed  in 

Equation  h ,16  and  ^,17i  the  following  ten  •  si  ould  b° 

-'Vied  to  each  pair  of  equations. 

a )  Heave 

i/L-V      •     WIWI      and         1/Z   Zzo3    .    <zo    -    1)2  (4,19) 

-Z»    +    m  -Z.    +   m 

w  w  , 

where  the  first  term  represents  the  fri c  :":  ona]  and  f  rm 
drag  and  the  other  term  is  the  second  term  in  th  ■ 
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polynomial  fit  of  the  Z  (?0,<£) 


<j>=0 


b)    roll 


-k.+I 

D        XX 


1/3!  f'^x3  +  W.BG)      (<j>--Si)3  and 

-v  +1 
fj      xx 


(4,20) 


l/3»  V/.?G 
-k.+I. 


"D 


XX 


t 


Where  the  first  term  is  again  the  contribution  of 
frictional  and  form  drag,  the  second  term  is  the 
contribution  from  the  polynomial  fit  of  the  curves  of 
stability  and  the  last  is  the  following  term  in  the 
expansi-on  of  GB  sin  ©. 

In  th  f3  case  3  *i'r^~>rr^l   th  q    vs  1  5  d  *  tv  of  th  s  6j    -  ■  *  ^>~i 
for  the  excitation  is  even  more  questionable  due  to 
the  small  orders  of  magnitude  involved  in  the  identifi- 
cation.  One  way  of  obtaining  a  better  approximation 
is  to  consider  the  effects  of  the  presence  of  the  wave 
ir  the  damping  terms  as  is  done  for  the  cose  ef  roll 
by  Tamya  [22]. 

Those  terms  would  then  look  like 

a )  heave 
\     .  (w-  jl      and   1/2  Z  a  ,  (w-  %)  •  \    ■   %   |        (4.21) 


-Z'+m 
w 


■Z.+m 

w 
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b)    roll 

\  (p-^}         and      l/2    kpa      .     (r~     ft'l1"     k\         (4*22V 


-k-+I  -k«+I 

p      xx  p      xx 


where  -\  -      ?*■ 

5t 


ft- 


6* 


j 


6t 

When    considering  the   non-linear   coup]         heave   and 
roll   mode,    all    the   coupling   terms   as     lentioned    in 
Chapter   3   bad    to    be   added    fro*"   which    only   tl  "■    pure 
hydrostatic    terms    :ould    be    pr    ~:  '  :ted    by  using,    for   example, 
the   method   presented    in    Ref   r    r.c     "r 6 ,    7  J.      ,T"~  ~    3£ 
objections   would         plj    wit]  pect   to   J"'   ;     >xci       bion. 

In   summary,    the    excitation   due    to    tl    =    :    :ti  3r.    of 
plane   progressive  waves    can   be   calculated    cjjiii   -      eci  rately 
using    i'v  s   Haskind   relati  )r.  ',    being    the   action    From   the 
ilif  fraction   potential    calcul;   bed    for   each       '  ttion 

two   dimension  a]    techni ,_  I    t)  e   resul  :  ar  t     '     'ce    s  id 

moment    obtair  3d    by    ini     "\"~  '  '    n    of   those    two    j  :     »jn  ■'       &J 
results    along  the    ler  -t]     -"r   the   vehicle.      Due    to 
computational     iifficul  '- :     ■  tl  method  d   which 

gives    results    :     1  '  tc    1       less     iccur    be    :  "  '         f  ;rs 

the   advc  ita  ;e    of    its    simplicity.       \not]   i  • 

-<-  5_ori  2(j    v  .fore    is    t]  tat.  this   later  method  ■       - 
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th  e  exciting  forces  and  moments  in  terms  of  instantaneous 
values  of  surface  displacement  and  wave  slope  which  are 
easy  to  obtain  during  the  tank  experiments  setup  to 
measure  displacements  and  velocities  that  will  "be  used 
in  the  identification  procedure. 
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CHAPTER  *> 


EXPERIMENTAL  PROCEDURE 

A  model  of  the  "mariner"  type  hull  form,  whoso 
characteristics  are  giver  in  Table  5»li  was  tested  in 
regular  waves  measuring  the  displacements  and  velocities 
in  heave,  roll  and  coupled  heave  and  roll,  and  also 
the  surface  elevation ,  v/ith  the  purpose  of  identifying 
the  hydrodynamic  coefficients  in  those  modes  of  motion. 

TABLE  5.1  -  MODEL  CHARACTERISTICS 

scale  (model/full  size)         l/96 

1  or.gth  (L)  5.50  ft, 

beam  (3)  .79  ft. 

draft  (T).  .28   ft. 

weight  (WQ)  ^6.05  lb. 

mass  (m)  l.'i-3  lb.r^c2/ft. 

The  model  was  placed  in  the  mid -length  of  the  tank 

with  its  centerline  perpendicular  to  the  centerline 

of  the  tank  such  that,  using;  the  coordinate  svstem  as 

defined  in  Chapter  2,  the  direction  of  propagation  of 

the  waves  generated  at  the  wavemaker  in  the  beginning 

of  the  tank,  coincide  with  the  negative  x   =  v  direction/ 

3 

The   model   was    placed    at   the   mid-length    of   the    tank 
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as  a  solution  of  compromise:  not  too  close  to  the 

wavemaker  in  order  to  avoid  the  notion  of  the  diffracted 

waves  from  the  model  upon  reflexion  on  the  wavemaker, 

not  too  far  from  the  wavemaker  in  order  to  avoid 

interference  from  the  walls  of  the  tank.   Reflexion  on 

the  opposite  end  of  the  wavemaker  was  of  much  less 

concern  due  to  the  existence  of  a  Hbeach"  composed  of 

stainless  steel  shavings  with  a  reflexion  coefficient 

of  about  .05 i  for  normal  wave  lengths. 

Tests  were  run  at  three  different  frequencies. 

First  it  was  decided  that  those  frequencies  would  he 

half,  twice  and  the  natural  undamped  frequency  for  b<^+1"' 

hea.v  3  and  roll,  where  the  natural  undamped  frequencv 

is  defined  as  for  a  linear  system 

CO  =J   k/m         where   k  =  -Z   -  -k<k 
n  zo     t 

m  -  -Z'+m  -  -K«+I 

w       P  xx 

For  roll  the  measured  natural  period  (from  free 

oscillations  of  the  model  in  calm  water)  was  T   =1,52 

n 

seconds  or  f   =  , 66  cos.   Therefore  tests  were  ru^  at 
n 

f1   =       .33    cos  A-    gT2/2iT  -   ^7.36    ft. 

f     -      .66  cos  =   1 1, 8k   ft. 

f      =1.32    cos  =      2.06    ft. 

3 
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For  heave  the  measured  natural  period  was  T   =  .78 
seconds  or  f^  =  1.28  ops.    For  the  tests  to  he  run  at 
twice  the  natural  frequency  f  =  2.48  cps,  it  would  be 
too  close  to  the  maximum  at  2,50  cps  permitted  by  the 
wave  maker.   Therefore,  it  was  decided  to  run  the  tests 
in  heave  at  one  fourth,  half  and  at  the  natural  frequency, 
or 

f.,     .32  cps  A  -  49.92  ft. 

64  cos  =  12.48  ft. 

=  3.12  ft. 

or  the  tests  in  coupled  heave  and  roll,  it  was 
;i(3.ed  to  run  at  the  natural  frequency  in.  !  ave  r^d 
at  half  and  at  the  nature.1  frequency  in  roll,   therefore, 
for  the  coupled  mode 

f,  =   .33  cps  X3  =  a 7 . 3 6  ft. 

f  =   .66  cps  A  =  11. '84  ft. 


1 

f  = 
2 


f  -  1.28  cps 
3 


f  =1.28  cps 


A  =   3.12  ft. 
3 


The  quantities  measured  in  each  test  wen 

a )  heave 

surface   ii^o1    ci  ment  ${t) 
heave   displacement      z0(t) 

heave  acceleration  z    (t) 

o 
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The  reason  for  measuring  acceleration  "nstead  of 
velocity  was  due  to  the  impossibility  of  obtaining  a 
translational  velocity  tr;  nsducer,  and  to  not  being 
'■hie  to  build  a  reliable  electronic  integrator  with  the 
material  available 

b)  roll 

surface  displacement   ~^(t) 

roll  displacement      G)(t) 

roll  velocity  p(t) 

A^rain  wave  slope  should  be  measured  in  this  mode 
ir.^tr"?d  of  surface  elevation  but  the  probe  used  could 
qv,  ]_y  rneasure  ){ t ) 

c)  coup]  2d  heave  and  roll 

surface  displacement   -^(t) 

v,  gave  d  is plac  ement    z  ( t ) 

h  e  a  v  e  a  c  c  =  3  2  r;  1 1  i  on    z  ( t ) 

o 

roll   displacement  &  (t) 

roll    velocity  p(t) 

The     >ame    observati    ns  as    ma  !~    :  n    a)    ai  lb);  vov^ 
p  l  v « 

i 
Description    of  probe  '  sur 

simulti      "    usly   surface     lisplacemenl    \         wav       ;]   >pe 

3    available    in    Reference    C9'''l» 
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Surface  displacement  was  mea.sured  with  resistance 
type  wave  probe  formed  by  two  parallel  sta.in]  ;  ;s  steel 
wires  of  l/8"  diameter.   Since  it  was  necessary  bo 

obtain' the  wave  displacement  at  the  vehicle,  the  ideal 
position  for  the  wave  probe  would  bo  along  the  centerline 
of  the  model.   Since  at  this  position  the  probe  would 
be  subject  to  interference  between  the  model  and  the 
tank  walls,  it  was  decided  to  put  the  probe  in  integer 
number  of  wave  lengths  ahead  of  +1  =  model.   Those 
distances  were  calculated  based  on  the  relation 
v  _  oj3  /g  or  J\=  gT  /->T  and  this  procedure  introduced 
+  v,a  first  error  on  the  measurements  since  F   r  Finite 
water  depty  :  the  relation  be    n  vave  nui  1  it   and 
frequency  is  C^3j 


<*>2  -  k  tanh(kh) 


(5.1) 


where  h  i c:  the  water  depth. 

For    sxample,    considerir      h    =   !i   ft, 
f   _  to  /2n  (^    -   ^h  )  //^~° 

.5    cps  .Qop 

.65    cps  .025 

7.Q5   ct  ~    0 

From   the   at  )ve    :  t    is    seen   that    for   the   ] 
frequencies    f<    .5   cps    tl  i^        ror   be      mes  '  ■  h  ble. 
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Heave  displacement  was  measured  with  a  linear  variable 
differential  transformer  (LVD?)  connected  to  the  heave  rod 
assembly.   The  heave  rod  slides  on  a  sleeve  provided 
with  ba] 1  bearings  to  reduce  friction.  Even  though,  the 
friction  was  found  to  be  quite  high  (comparing  with 
experiments  had  with  a  heave  rod  assembly  provided 
with  air  bearings)  and  therefore  there  will  be  solce 
error  involved  in  the  identification  of  the  heave 
d  a  mp  i  n ,;"  pa  ra  m  e  t  er s  . 

Heave  acceleration  was  measured  with  an  ^ccelero- 
meter  placed  in  the  model,  as  close  as  possible  to  the 
axis  of  rotation  n"n  order  J~o  minimize  in'  srference  of 
roll  i  n  the  me?  surement  of  heave  ace  il   ?r;  t:  in. 

Roll  displacement  was  measured  by  a  rotational 
variable  differential  transformer  attached  to  a  roll 
bearing. 

Roll  velocitv  was  measured  bv  means  of  a.  rate 
gyro  t  bo^ro'.'.'ed  from  the  CI  irles  Draper  Laboratory,  which 
was  placed  in  the  model  with  its  thr  3  axes  of  *oi  '" 
parallel  to  the  coordinate  -  xes  located  on    1}   =   jdel. 
This  rate  cyro,  used  in  the  Deep  Sub]      :e  Syste  is 
Project,  was  able  of  m  a:  irii     i  ular  rat  s  "   roll, 
pitch  and  yaw  simultc   i  u  ;ly,  up  to  rs  fces  of  ^0  dc     •/ 
■  ond  . 
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All  the  transducers  mentioned  above,  with  the 

exception  of  the  rate  gyro,  are  excited  by  an  A.C.  voltage 
of  5  volts  at  a  frequency  of  2*4-00  cps  provided  by 
carrier  preamplifiers  Sanborn  model  pQ05A  which  also 
received  the  electrical  response  from  the  :  -ansducers, 
generating  a  D.C.  signal  which  may  be  fed  to  a  chart  or 
tape  recorder  or  both.   The  amplifiers  also  provide 
n  one  of  amplifying  or  attenuating  the  electrical  signal 
according  to  its  level  since  both  recorders  work  between 

±2.5  volts. 

Since   the    excitation  necessary   for   the   rate   gyro 
•  .    26  V  A.C.    (f  =   ^00   cps)  ,    the   Ss nl  c  ■  ■     pr  lif i  ers 

could   not    be    used,      Inste    d    the    following   circuit  was 
use(3    (see    pig,    5.1),    since   a   D.C.    power   supply  was 
available    up    to"   three    amps: 


r 

24  V 
DC 


/  3 

400  Hz 
STATIC 

INVEP-TER 
2  4 


__  C  ^.47 

AF 


115  V 

\400Hz 


26V 

\400  Hz 


Pn       mr»D        ^    .    1 


The   26   V   'K)0   hz   excitation      ■  s    fed    to    bh      rate 
gyro.      The   output    from  the   rate    gyro    (!'-00  hz.    ph;    ?e 
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reVersing)  was  then  fed  to  a  demodulator  (PRD  type  808) 
to  obtain  the  D.C.  signal  proportional  to  the  angular 

velocity.  This  signal  was  attenuated  in  a  Sanborn  D.C. 
amplifier  before  going  into  the  chart  and  tape  recorder 
(see  Fig.  5.2) . 


FROM 
GYRO 


1 
115V 

^400  Hz 


3  4 

2  DEMOD. 

PDR  80  Q 


T 
X 


5/-F 


D]CT0     DC 

J 


AMR 


Figure  5.2 

For  convenience  of  future  calculation,  it  was  decided 
to  out  the  center  of  gravity  at  the  design  waterline 
(p.l]  stability  calculatioi  i    n  made  for  en  assumed 
KG   .28  -— ■  •  ~  center  o^  ;ravity  at  the  wa1  ■■'''■     :)  .   This 
was  performed  by  trial  ar  I  error  by  creatin   e  kn<  ■ 
moment  (with  the  model  in  the  upright  positioi  ,  a  g] 
weight  is  moved  transversally  a  certain  giver  ;:-'-:  ;e) 
an(^  r^easurin  ■•  the  resulting  inclination  which  '/as      red 
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with  the  value  calculated  in  the  curves  of  stability. 
'»fgj  s;hts  were  then  moved  "•  n  the  vertical  and  the  th^Is 
verified   again  .until   reaching  the    value    s       ■*  '       the 

curves . 

Calibration   o^  all    instrumentati<  '      in    tne 

tank   except   for  the   rate   gyro  which    was   performed    a1     fhe 
Draper  Laboratory. 

The   LVD?  was    first    zeroed    by   slidii         he   core    irrto 
the    casing  until   no    output   was    s]      t    an  the    chart   recorder. 
This    point    ,vas    mark  5d    on    the    sha  Ft    of    t}  =         r      •  -  --  '-^ 
■'.'"'  I'.h   niG  r^s   s  os  c  ed    3  "    <^n    6  2  ch    s  ]  d  ,_  .      Th  en    f  or  a.  ^c  t 

of   the   core   of   1"   and    for  a    -^  ven   s  t  ;        s  ti  on    o  f    I '        ^re  — 
a  hit)]  if  is  r.i  >    bh      oe;     ief  1     itior      i     bhe    '■ '  "       )vd    c  was 

noted    (1"   =    10'  mm   -    5^    "v/mrn   attenuate  i     f  ')  ,      For 

the  wave   probe   and    rots  bional   di fferentia]    tr    isfo 
the   orocedure   was    similar.      For   the   rots    !       i!    tr;  nsfon  er 

brat i or    of   th*    heave       'n         I   :     11    bearing  was    checked    later 
with   the    instrumentation  abo:   ?& . 

For   t1:       ;ali1        bic  '  the   acceler  ■     :    : 

checked   v.' it]     it   ii      bhe  1        '  '       bal . 

was    th -"-     I     rn  1 3    90     '  nd    ii       lu.1]        :alo  valu 

(1,5     ■  <>  .'■■  8  ,  3     rt .  7"  20  .  P  )       '  "ull       :ale  6     '-■■■■■ 
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th  e  pen  (25  mr" )  .   Then  the  zero  suppression  control  on  ' 
preamplifier  was  adjusted  to  bring  the  output  to  zero. 
For  the  rate  gyro,  the  whole  assembly  (rate  gyrof 
inverter,  transformer  and  demodulator)  was  placed  in  a 
rotating  table  to  which  a  known  constant  angular  velocity 
was  given.   The  output  voltage  was  then  measured  with  a 
high  precision  D.C.  voltmeter.   It  was  noted  that  'the 
presence  of  demodulator  introduced  some  hysteresis  on 
the  output  as  well  as  some  non-linearity  above  an  angular 
velocity  of  20  degrees/second.   On  the  average  the  output 
signal  was  of  , 136v/degree/second  D.C,  while  for 
the  rate  syro  alone  5  ;:  was  . l^-7V/degree/second  rms . 
Afterwards  in  the  tank  a  constant  an  ular  velocity  was 
simulated  by  means  of  a  potentiometer,  the  output  voltage 
r^ad  with  an  oscilloscope  and  the  pen  deflection  on  the 
chart  recorder  noted. 

The  results  of  the  final  calibration  for  all  the 
ii  struments  were  as  follows: 

•"eve  probe  1"  =  10  mm  (atten,  factor-  =  10) 

heave  rod  1"  =  10  mm  (atten.  factor  =  5°) 

accelerometer     1.5  g  r~  25  mm  (;  bten.  factor  =  50) 
roll  bearing        10°  =  10  mm  (atten.  factor  =  200) 
rate  gyro      ?0°/sec,  =  10  mm  (atten.  factor  =  200) 
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Froin  the  above  so^e  other  mistakes  car.  be  noted. 

First,  full  scale  (25  mm)  corresponds  to  2.5  volts.   For 

nodel  of  3  •  375  inches  of  draft,  heave  displacements  of 

more  than  1  inch  (  1/3  T)  can  hardly  be  expected.   Therefore, 

the  output  signal  on  tspe  will  vary  bet1     £  1  volt  at 

most.   For  the  acceleration  if  it  is  assumed  that,  as  for 

linear  svstems  z  —   lc^z      the  amplitude  of  tht*3  3c^olpr^ tion 

a      a       * 

woul 1  be  for  the  highest  frequency  (f  =  3  .9P  ops 

-  Q .  0^-  sec ,   )  z   -  64  .  64  inch/s ec .  3  or  5  .  3P  f  t ,  /s ee  ;  s 
since  full  scale  was  1.5  g  =  '^  •  3  "p".  ^s^c.   the  signal 
or  tape  would  vary  between  ±  .25  volts.   For  the  lower 
pr'oc! uenc ]  es  th  2  accel  o""n  !_^  ^"^  ou  tnut  was  ,',;'  cdl  p  i  srceiDt^b"'  ° 
qt  , . '  ■  >  • :-  .  ••       •  not  jsabl e  ■f'or  th e  d  i  ^it  i  z;  '      s in c ^ 
output  signal  approached  the  ^oisc  level. 

Similarly  for  roll  displacement  and  velocity  and  for 
the  wave  probe,  the  signal  level  was  quite  low. 

Before  running  the  tests  calibration  ";  nals  were 
"'-■>!+  into  the  tare  recorder  in  order  to  be  ab!  rA  to  sca.l  g 
1 -1,-.  g  .i •-! -j-g^  -\- q  r-. ^■pT'ec t  v3  1  u  ss  o  F  lis ol  acements    ;  vel  cities 
once  the  digitization  being;  performed^  About  50  ?i  .  of 
taoe  were  recorded  wi  th  l~  o  >e  ca]  ibratior  si  it  lis  at  p 
->-> g  >->  o r d  e r  s n e e  i  o  f  7i  i  n  o  ^  /s  <"■ c 

Then  the  nine  tests  as  rr   tioi  vl  before      run 
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recording  about  50  cycles  for  each  run.   I>  between  runs 
enough  time  would  elapse  in  order  for  the  :rodel  and  t\  \ 
water  to  calm  down  (zero  initial  cor.']  it  ion."  )  .   On  tape 
the  interval  between  runs  was  about  20  ft. 

The  waves  were  generated  by  a  wave      c   constituted 
by  a  hinged  paddle  actuated,  by  an  electro-hydraulic 
power  system,  which,  was  controlled  by  electrical 
sinosoidal  signals  generated  by  a  low  frequency  wave 
oscillator,  through  a  hydraulic  servo  valve. 

The  data  was  digitized  from  the  me.  •   bic  tape  into 
computer  cards  by  means  of  the  SAI  680  ar  '  :   c   puter, 
EAI  693  :'  ts  conv  jrsion  system  an       L130  CEM   '  Lta] 
computer  of  the  Depari    :  of  [V!ecl  <   nical  Sngineeri]  ;, 
During  ligitization  some  more  defic]   lies  were- noted , 
besides  the  low  level  of  the  signals  en  tape  (the  FAT  680 
works  within  a  voltage  level  of  ±  10.5  ',r,~ 1  ''  '    '■'"    nor  :al 
operation;  when  digitized  thosf   •'    Is  to  ■' 

standard  precision  real  valus  of  '-    1,05) •   First  tl  3 

;netic     :  '■  e  :  ,  on  the  po]  ::         recc     ,  a  D.C, 
signal  e°  about  5«^  voli       i   rigi   ;<     '  i 

i         ' ;  p  y  •■  ■  i  ■-  •    i(j  fi       ■  •  '      1  •  1  g  1        q  '      '      -       ■■,,.•'       -      -j    ^      C  OU  1  d       *    '  '  |  70 

the   heads    of    tl   2    tape    c   c  ) :  3  ir  w]   jr    on    "  »]  s.y" ,    thus 
introducing   errors    on   the   values    of    tl  ;ur  ;  ". 
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Also  the  runs  were  too  short  to  ho  possible  to  net 
up  the  values  of  the  potentiometers  conveniently.   These 
potentiometers  were  used  to  take  off  any  bias  tha.t 
existed  in  the  measurements. 

For  purposes  of  digitization,  these  circuits  were 
built  in  the  EAT  680  analog  computer  (refer  to  Appendix 
k   and  Reference  r?5~l  ^or   details  in  the  circuitry).: 

1 .  -  Control  of  Sampling  Pete 

A  Wavetek  116  VCG  oscillator  is  used  to  generate 
a  square  wave  of  frequency  equal  to  the  sampling  rate. 
This  signal  is  obtained  in  analog  trunks  0  or  3-   From 
■"v:c  trunks  the  signal  goes  into  the  input  of  a  -■         p-  vir 
in  the  interface  tray  of  the  dialog  rows.   The  :  rtput, 
which  is  obtained  in  the  same  interface  tray  in  the  logic 
row  is  fed  into  a  differentiator  to  obtain  a  spike*   The 
output  of  the  differentiator  is  in  turn  fed  into  the 
analog  to  digital  converter  (ADC)  sample  control 
terminals  on  the  first  strip,  preferably  on  '    '  Mal] ", 
This  signal  may  also  be  fed  to  a  BDC  count  r  _:  \     C  lu  n  )1 
to  check  the  operation  of  the  circuit, 

-  •  ~  S  i  rn  a  1  F  r  o  c  ^  s  s  i  n  g . 

From  the  tape  recorder  the  signals  are  ;:*  ?d 
into  analo   trunks  1^  to  20,  located .below  the  Wavet 
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oscillator,   From  the  trunks  '}    ise  output  is  situated  in 

second  r^^r:"  "'",    ■•    ■     .  •  <.,  .     -   r  together 
with  the  output  of  a  potentiometer  (to  take  off  any  1  ias 
on  the  si  jnal)  .   The  output  of  the  summer  is  Fed  into  s 

ntiometer  "~v"1<  frorr  "■   potent ic   "    '  '   -   amplifi  c 
(input  marked  with  a  10).   r"-  e  purpose  of  the    '  r.tio- 

" ter  is  to  attenuate ,  if  neces  i   ~,  ' '   ■ "  •  ■  ]      "  :  ch 
ir  amplified  ten  times  in  the  amplifier.   The  signa] 
i  rrespondent  to  the  measured  ace  ]  *;  bion,  "•  stead  of. 
go  ins  into  an  amp]  5  fier.  is  fed  into  an  '  _  :   rator  to 
convert  '  '"  to  velocity,   The  ou1   t  of  fc'      ]  '  "'er 
^-.  *  v-.  -!-  2 r"  'a  tor  i  s  tin  er  fed  int  c  th<  '  DC  tr      ( 1  oca.  ted 
■  -.-  bh e  s  ec ond   brin)  '   I    '  '  '   d        '  ■    ~   ■ 
+  ■  ry  -,  ■  ■ -.  l,-  the  5  is  d]  3  ,r  trunks  for  d  is  q]  av  5  '  '   the  ^FS 


DU390T  oscilloscope. 

•  -  - ' -'  •-  - J  --  -  -  -, 

From    t  ;rmii  '^ n     ';    for   >  pie ,    in  the 

3  0]  tro]    timer    (in   column   02   of  the   ]       :        ■      )   a  wire 

connected   to   the    "scope"   trunl 

TRC 

•       |  ] , ,-    p-,.,,-1-     ;-tr"ip    in    order   to    ol  '     :_  lisplay   5      any 

ana.]  o u   \   )d  e    c  o t  t  ~ '  ■ i  n         \    oil     is    i     I       iuc ed    in    ths    " v"  i    " 
t         '       ]    3     :ated   close  terminal     3  '  i"bove. 
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Digitization   war-    initiated     '  IB!     1130   by  a 

computer   program   built   according   to    the   hyl  rid       ser's    5"ide 
du bl  i sin  ad    by   t b °   D ere  r tmer "t    0  f   Me 0 '  ~  •  ica!!    E  -.->•■  *  - 

using   existing  routines.      This    c ber    pr  is 

presented    in  Appendix   h  . 

Tin o    outDut    from   the  r^  *  "^  t]  r"  ,_ "         was    t  ; '  2T  '     '  0 

c  0 ■•  01  iter  e 2  r d  s  . 

The   selected    sampling  ra^  ■    '   ■  nty   s;    -  ■'  >s/ 

cyole    for   fifty  cycles    for   ^i"1    runs,       '    ch    n    freou,  ?nt 
sample   will   allow  a   method    of   linear    Lnt        ol<  '  ion    tt    be 
i^'^j    in   the    identification   of    th<  1    ~r"s . 

t -v-i    ■  •  1  !...■<•  yi •  -   th '     '     ]  '     Doinijs    '  '        n      p ■  t ri 1 

procedure   were    ou  bl.ined ,         I  1st;      g  mad 9 ,     '  '       y 

to    inexperience    in   this   kind    of  'im      bs ,     \rere   a] 

pointed    out.      If   time    permil  bed     t)  \     should    be   repea  in 

■  •■■  i  3r   to    obtain    bett<  ?   data. 

Digitization   w   s         plained    in  s       ■  '  ler 

4- q    r=  "i  1  ow   ^    s om^wha-!       '*  '  '  e  c      -      •  1  "  )*p     (g  J  ■      ^  ■  ■     m  ■  '   ■ 
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CHAPTER  6 

IDENTIFICATION  OF  SHIP  MOTION  MODEL  PARAMETERS 

Most  of  the  available  techniques  of  parameter  identi- 
fication have  been  developed  for  linear  dynamic  systems 
[l].   Some  of  those  techniques,  assuming  the  linearized 
equations  of  motion  to  be  a  valid  description  of  the 
behaviour  of  an  ocean  vehicle,  have  already  been  applied 
in  the  determination  of  ship  motion  parameters  in  heave 
and  pitch  namely,  the  frequency  response  Q27I  and  the 
step  response  methods  C28~J. 

In  the  case  of  a  non-linear  model,  two  approaches 
are  open,  both  yielding  approximate  results. 

The  first,  making  no  assumptions  about  the  character- 
istic of  the  uncertainty,  assumes  a  mathematical  model 
for  the  system  and  then  comparing  the  output  of  both 
model  and  systems  to  the  same  input,  searches  for  the 
structure  of  the  model  that  minimizes  a  function  of  the 
error  betv/een  the  two  outputs.   An  example  of  the  fore- 
going is  the  model  reference  technique,.  For  the  system: 

y  (t)   =  f  (u(t,t)  ,t) 
s        s    o 

y  (t)   =  f  (u(t  ,t)  ,t)  (6.1) 

m       m    o 

V      =  F(y  -x   ) 
s   m 
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where  u  vector  input  to  both  system  and  model 
v   -  vector  output  of  the  system 

y„  -  vector  output  of  the  model 

^m  c 

fg  -  vector  structure  of  the  system 

f   -  vector  structure  of  the  model 
— m 

V   -  scalar  function  to  be  minimized. 

A  simplified  version  of  the  above  will  be 
y_s(t)  =  f(u(to,t),  t,  j^) 

zm(-t)  =  fiuCtgtt),  t,  p_m) 


(6.2) 


where  the  values  of  p_  are  unknown  and  the  values  of  o 

will  be  determined  by  minimizing  V,  which  is  often  of 

the  form 

v(£m,tQ,t)  =  _i_  ft  [ia(t)-yto(u);]C2£e(tf)-2m(x)f  dt       (6.3) 

o   o 


and  for  optimum  p_   (best  approximation)  6V  =  0.  The 

m  6o 

J~m 

vector  p_  thus  found  will  only  be  optimum  in  the  interval 

t  <  V  <t   and  may  not  apply  outside  that  interval.   Also 

it  may  not  be  true  for  non-linear  systems  that 

lira   2=2  (6.^) 

t— •>«>   m    s 

as  one  might  expect  from  linear  systems.   The  optimum 

value  of  p_   in  Equation  (6.3)  could  be  found  by  a  maximum 
m 

gradient  seeking  algorithm  (first  order  gradient  method) 
or  by  a  second  order  gradient  method  £29] 
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The  second  approach  is  to  use  linear  estimation 
techniques,  such  as  a  Kalman  filter  to  the  linearized 
form  of  a  modified  state  representation  of  the  system. 

Considering  a  linear  system  driven  by  noise,  Bryson 
and  Ko  [29,  Chp.  11  and  12]  present  the  following  expressions 
for  the  estimate  of  the  system: 

1)  Linear  Multistage  Systems 

x    =  F.x.  +  G.w 
-i+1   ~i_1   ~i~i 


x   -   state   vector  E[x     '   =  x 


E[x   -x   ]   [x   -x   ]     =  M  (6.5) 

~0   — o        ~~o   ~~o  o 


w   -   environment  noise  E|_w]   =  w 


T 


E[w. -w. ]    [w . -w  . ] L   =  Q . 6 .  . 

_  _     T 

Erw.-w.l    [x    -x    ]      =    0 
u    1  — iJ    ' — o  —o 

F.  G.  -  transition  matrices 
—l  —i 

and   z.  =  H.  x  +  v. 
~i   -i  -i   -i 

z_  -   measurement  vector 
H-  -  measurement  matrix  (6.6) 

V.  -  measurement  noise    E[v. ~]  =  0 

E!~x  -x  1  v  =0 
— o  ~0  ~ 
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Assurning  the  sequences   i  of  the  different  vectors 

to  represent  Gauss-Markov  random  sequences,   they   present 

the  solution  for  the  estimate  of  the  state  x  -  x  asi 

x.  =  x.  +  K. (z.  -  H.x. ) 
— 1   — 1   —l—i   — i— l' 

ii+1   =  Hi  *i  +  Gi  Wi 

K.  =  P^i^i"1  (6-7) 

r   -1       T   -1   -,-1 

Ei  =  CMj.   +  H.  Si  Si]   = 

Mi  -  MiHiT(H.M.H1T  +  R.)"1  H.Mi 

T         T 
M.  ,  =  F.P.F.   +  G-Q.G. 
—i+l   —l—i—i    —l—i—i 

where  P.  is  the  error  covariance  matrix  of  x.  after  the 
—i  -i 

measurement  while  M.  is  the  error  covariance  matrix 
before  the  measurement.   Similarly  x.  is  the  estimate  of 
the  system  after  the  measurement  while  x.  is  the  estimate 
of  the  system  before  measurement. 

2 )  Linear  Continuous  Systems  with  Continuous 
Measurements 

The  mathematical  expression  for  this  type  of 
system  is  obtained  by  letting  the  time  between  stages  go 
to  zero.   Therefore,  for  the  system 


x  =  Fx  +  Gw 

(6.8) 

z  =  Hx  +  v 
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with   F  =  lim   F(t-  +  a  ,  t- )  -  I 
A-0 7 


(6.8) 


G  =  lim  Gi 
a— 0  — 


Under  the  assumption  that  x  represents  a   Gauss-Markov 
randon  process,  the  expression  for  the  Kalman  filter 
will  be 

x  =  Fx  +  Gw  +  K(z  -  Hx) 

K  =  P  H   R  L 

£  =  !  p  +  £  £T  +  2  i  f  -£HT  BT1   h  p 

p  <*o>  =  Po 

If,  for  the  system  as  specified  in  (1.11)  to  (1.13) 

k*  =  f'CxIutp'J  +  *'  d.ii) 

p_  =  0  (1.12) 

z  =  hV  +  v  (1.13) 

the  above  technique  is  to  be  used,  it  is  necessary  by 


first  to  form  the  extended  state  vector 

X  = 


The  new  state  equations  are  as  follows: 
x  =  f(x,u)  +  w   S[x(t0)]  =  xQ 
z  =  H  x  +  v 
where  in  addition  to  6.10 


(6.10) 


(6.11) 


w  = 


w 


0 

u —  J 


(6.12) 
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f  =   f  •' 

.0  j  (6.12) 

H  =  [H'   o] 

Now  the  state  equations  are  linearized  in  order  to 
use  Equation  (6.9)  or: 


x  =  f(x,u)  +  w  +  K(z-H  x)5x(t0)  =  xQ 

T  -1 

K  =  P  HXR 

P  =   6f  P  +  P  6f  T  +  Q  -  P  HTR"1HP;  P(tQ)  =  P 


(6.13) 


6x        6x 
The  next  simplification  will  be,  considering  that  in 
practice  most  measurements  are  discrete  in  time,  to 
approximate  the  above  differential  equations  by  a  set  of 
difference  equations  and  use  the  formulation  for  a  linear 
multistage  processes  or, in  other  words,  for  the  system 
below  : 

-i+1  =  -i  +^*£(£i>iii)  +  A'^i 


(6.1*0 


z .  =  H  x.  +  v 


The  resulting  filter  will  be,  assuming  for  convenience 
that 

eIa'w^]  =  0 
L+l  =  ^i  +  A'£(Si.lii)  +  *'Ki+(!i  -  H  lL)  (6.15) 


A 

X 


T       -1 
K.  =  P.H  (R  .A) 
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P.  =  Hi  -  A'HiHT(H  Mi  H?A+  R-  A  )-1  H  M. 

N  m  (6.15) 

M.  _  =  _6 (x.+A»f(x,  ,u,  ))  P.  6_(x.+A-f(x-  ,u,))  l+   Q-A 

T+l       T 1       1   1  1  "1 1       X   1  

x  x    6x.  6x. 

-l  —l 

where  R  and  Q  are  constant  matrices.   If  the  value  of  Z\ 
is  sufficiently  small,  this  formulation  will  represent 
closely  Equations  (6.13)  ai"td  will  be  much  easier  compu- 
tationally. 

In  the  above  formulations,  the  partial  derivatives' 
may  be  evaluated  in  a  nominal  path  or  if  greater  accuracy 

A 

is  desired  they  are  evaluated  at  x  =  x. 

Bryson  and  Ho  [29]  mention  several  authors  who  have 
achieved  success  using  Equations  (6.13)and  (6. 15)  for 
non-linear  systems  but  at  the  same  time  make  the  warning 
that  convergence  may  not  be  obtained  if  the  initial  guess 
is  poor  or  if  the  disturbances  are  so  large  that  the 
linearization  is  inadequate  to  describe  the  system. 

From  the  foregoing  Equations  (6,2)  and  (6.3)  , 
Equations  (6.13)  and  Equations  (6. 15)  it  was  decided  to 
use  Equations  (6,15)  clue  mainly  to  computational 
convenience.   Due  also  to  the  poorness  of  the  available 
data,  as  mentioned  in  Chapter  5?  it  was  decided  to  abandon 
the  identification  in  heave  and  coupled  heave  and  roll 
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and  concentrate  exclusively  on  the  identification  of  roll 
motion  parameters  whose  data  at  a  frequency  of  excitation 
of  .66   cycles/second  (roll  natural  undamped  frequency) 
seemed  reasonably  good. 

As  mentioned  in  Chapter  U-,    Equation  (Jj-,1?)  the  state 
equations  for  roll  motion,  ignoring  for  simplicity  the 
dependency  of  the  system  on  the  derivative  of  the  input 
are: 

P  =    Kp   p  +  [  K4    -   W.BG  ]  (<j>-?y)+ 

-K'  +  I     -K'+Ivv 
p  xx    p  xx 


-K'  +  I 
p     xx 

W .  BG 

"Ki+Ixx 

+  [1/2 !Kp    ]  p|p|  +  [1/3 !K  3  + 


-K'  +  I 
p  xx 


(fc.l?) 


•K«  +  I 
p  xx 


1/38W.BG  1  (4-fv)3  +  -ri/3!W.BG  1  (b* 

-K'+I  Jy    .      -KA+I       T- 

P   xx  p   XX 


<f>- 


Non-dimensionalizing  the  above  equations  with  the 
frequency  of  excitation   the  following  is  obtained: 


t> 


CJ 


-K'+I 
p   xx 


1 


+ 


K 


-K'  +  I 
p  xx 


1 

0? 


W .  BG 


-K'  +  I 
p  xx 


OJ' 


*"fyl 


+ 


W .  BG 


-K'  +  I    CO 
Kp  xxx 


^ 


+ 


-K'+I, 


xx 


(6.16) 
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1/3 !K  3   1_  + 

-K'  +  I  -  w2 
_  p   xx 


1/3IW.BG   1_ 

-K'  +  I      co' 
p  xx 


(<HJ 


l/3'.W.BG   1 


L"Kp+Ixx 


u) 


I 


6U 


U>\ 


Based  on  results  of  Appendices  1,  2  and  3,  the 
following  non-dimensional  parameters  will  be, for 
uj  =  ^-.1^5  sec."1: 

=   .01  x(3) 


xW 


K 
o 

1      = 

CO 

-.00^2 

"Kp+Ixx 

K 
P 

1    = 

z? 

-.989 

-K'+I 
p      xx 

W.BG 

-K'  +  I 
p      xx 

1    = 

U>2 

2 . 7*j4 

1/2  !S  p 

= 

-.126 

-KP+Ixx 

1/3  IK    3 

1    = 

3.123 

-K'+I    ' 
p       XX 

— 2" 

to 

=  10  x(5) 


x(6) 


(6.16) 


(6.17) 


Also  and  to  have  the  values  of  all  the  state  variables 
between  ±1,  it  was  chosen  that 


1 

p_    = 

CO3 

.1  yd) 

.1   x(l) 

4    = 

to 

.1  y(2) 

(6.18) 
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=.lx(2) 

(6.18) 
=.lu(l) 

The  complete  set  will  then  have  the  form  ; 
.ly(l)  -  .01x(3).lx(l)  f  [x(4)-10x(5)]  [.lx(2)-.lu(l)]  + 
+  lOx ( 5 ) . Ix ( 2 )   + 

+  x(6).lx(l).l|x(l)|  -  [l0x(7)-.167  10x(5)] 
[.lx(2)-.lu(D]  3~  .167  10x(5)[;ix.(2)]  3 
.ly(2)  =  .lx(l) 
or,  rearranging 

y(l)  =  .01  x(3)  x(l)  *  \k(4)    -  10  x(5)]  [x(2)  -  u(l)]  + 
+  10  x(5)  x(2)   +  (6,19) 

+  .1  x(6)  x(l)  ABS  x(l)+[.l  x(7)  -  .0167  x(5)] 
[x(2)  -  u(l)]3   -   .0167  x(5)  x(2)3 
y(2)  =  x(l) 

Before  using  the  data,  this  model  v/as  tested  by  si- 
mulating noisy  measurements   —  only  x(l)  and  x(2)  were 
measurable  in  the  simulation  —  and  the  control  u(l), 
with  a  format  similar  to  the  one  obtained  in  the  data 
collected. 

Also  simplified  forms  of  the  state  equations  were 
tested.  Their  expression  were  : 

No n~ linear  equations  of  motion  to  Identify  linear  terms 
y(l)  =  .01  x(3)  x(l)  +  [x(4)  r  10  x(5)][x(2)  ~  u(l)]  + 

+  lo  x(5)  x(2)  +  .0126  x(l)  AES(  x(l)  )   + 

*     (6.20) 
+  [.03123  -  .0167  x(5)][x(2)  -  u(l)]   + 


+  .0167  x(5)  x(2)3 
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y(2)  -  x(l)  (6.20) 

Linear  Equations  of  Motion 

y(l)  =  .01  x(3)  x(l)  +  [x(4)  -10  x(5)]  L>(2)- 

u(l)]  +  10  x(5)  x(2)  (6.21) 

y(2)  =  x(l) 

Equations  (6.19),  (6.20)  and  (6.21)  were  tested 
varying  the  following  parameters: 

1.  Initial  State  Estimate  -  The  "true  values"  of 

the  parameters  (x(3)  through  x(7))were  set  in  the  simulation 
at  the  values  given  by  Equation  (6.17).   The  initial 
estimate  was  set  either  zero  or  Equation  (6.17). 

2.  Initial  State  Estimate  Error  Covariance _P0  -  This 
is  a  diagonal  matrix  of  equal  elements  which  were  set  to 

0  ,  2.  and  5  . 

3.  Order  of  Magnitude  of  the  Elements  of  R  Relative 
to  the  Elements  of  Q  -  Both  these  matrices  are  diagonal  and 
of  equal  elements.   The  ratio  of  the  elements  of  R  to  the 
elements  of  Q  was  set  at  1 ,  2,  5  and  10. 

U- .      Noise  Level  in  Environmental  Noise  -  w  and  In 
Measurement  Noise  -  v  .  Set- at  0  and  1 

Then  the  collected  data  was  used  in  the  identification 
procedure  introducing  small  variations  in  the  parameters 
mentioned  in  1,  2  and  3  above  around . a  reasonable  set 
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obtained  in  the  simulation  to  obtain  the  best  possible 
convergence  and  also  to  test  the  effect  of  initial 
conditions  -  or  in  other  words,  the  necessary  accuracy 
in  the  theoretical  prediction  of  the  values  of  the  para- 
meters, to  obtain  convergence.   Some  of  the  results  from 
the  identification  are  presented  in  Appendix  5  while  the 
discussion  of  the  results  is  presented  in  Chapter  7. 

Summarizing  in  this  chapter  some  of  the  available 
techniques  for  parameter  identification  were  discussed* 
Attention  was  focused  in  those  that  permit  application 
to  non-linear  dynamic  systems.   From  those,  one  was  chosen 
that  having  perhaps  less  accuracy  in  the  results, is  easy 
to  implement.   The  equations  of  motion  for  roll  went 
then  through  the  necessary  modifications  to  be  in  agreement 
with  the  format  of  such  formulation.   The  structure  of 
the  model  and  the  formulation  were  tested  by  simulating 
the  excitation  and  the  measurements,  and  varying  selected 
parameters  to  obtain  a  better  insight  of  what  to  expect 
from  the  identification. 
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CHAPTER  7 

DISCUSSION  OF  RESULTS 
1*  Simulation 

The  first  equation  to  be  tested  was  Equation 
6,19.  The  noise  level  was  set  to  zero,  the  ratio  of  ele- 
ments of  R  to  elements  of  Q,  from  now  on  expressed  formally 
as  r/q  by  similarity  with  the  scalar  case,  equal  to  1, 
the  initial  state  estimate  error  covariance  matrix  P^  at 
[2]  and  the  initial  state  estimate  x0  at  zero  ( the  "true 
value"  wa3  set  in  all  simulations  at  the  values  expressed 
in  equation  (6.17)  ).  The  values  of  the  estimate  bj.ew  up 
and  after  200  measurements  the  computer  indicated  overflow. 
This  type  of  behaviour  with  no  noise  decided  the  use  of 
Equations  (6.20)  or  (6.21).  These  two  models  were  both 
tested  at  no  noise •  Equations  (6.20)  were  first  tested 
with  x(0)  -  Eq.(6.17),  r/q  =  1  and  P 0  =     [5]  .  The  error  on 
the  estimate  (XHAT  -  XX)  where  XHAT  is  the  estimate  and 
XX  is  the  "true  value"  of  the  state , remained  at  zero  during 
the  whole  experiment.  However  P(3,3)  (variance  of  x(3)  - 
linear  damping  coefficient)  kept  increasing  with  small 
oscillations*  The  other  elements  of  P  decreased  quickly 
to  less  than  01  •  The  same  behaviour  was  observed  for 
Equations  (6.21),  showing  that  the  influence  of  the 
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non-linear   terms   is    small,    at   least   for    the    no   noise 
condition. 

The    last  experiment   in  the   no  noise    series   "was  con- 
ducted with  Equations    (6.21)    Pq  =    [2]  ,    r/q   -    1    and 
(XHAT)0=    [OJ. 

The    error    (XHAT  -  XX)    for   all  variables   converged 
rapidly  to   zero  -    in  less   than  100  measurements   with  the 
exception  of    (XHAT(3)    -  XX(3))    whose   value    remained 
constant   throughout.   Again  the    term  P(3,3)    steadily 
increased;   while    all   the    other  variances   came    to  values, 

after  500  measurements,    less   than  .060    (initial  value  was 
2.0)    the    final  value    of  P(3,3)    was  2.202.    The    foregoing 
suggests   that  when  no   noise   i3   present   the    system  converges 
independently  of  initial  conditions,   with  the   exception 
of  Equations    (6.19)    and  with  the   exception  of   the    linear 
damping  coefficient  x(3)    whose   identification  seems  very 
difficult   to   obtain. 

The    rest  of  the   experiments  were   performed  with  noise 
level  equal   to    1  which   is  believed   to   be    quite    a  high 
level  due    to   the    scattering  observed   in  the   measurement 
vector  z. 

The    first    simulations  were    taken  with  r/q  =    1 
XHATp -    [0]  and   P^=     £3]    and    =   [5].    Both   simulations   blew 
up   at   about   the    140th  measurement.    At   that   time    p(3,3), 
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P(4,4)  and  P(5,5)  were  quite  large,  but  different  from 
one  case  to  the  other.  The  main  point  to  be  observed  here 
is  the  conjugate  effect  of  noise  and  not  accurate  initial 
estimate  . 

The  next  experiment  put  together  Equations  (6,20)  and 
(6.21).  For  both  systems  r/q  =  1,  P^r  [o]  and  (XHAT)Q  = 
(XX)  •  Since  initially  the  error  covariance  matrix  was 
set  to  zero,  it  increased  and  then  oscillated  about 
values  equivalent  to  the  ones  found  in  the  no  noise 
condition  (less  than  .050),  with  the  exception  of  P(3,3) 
which  steadily  increased  to  a  final  value  at  the  500th 
measurement  of  .249  in  both  cases,  with  tendency  to 
increase.  The  errors  on  the  variables  oscillates  around 
zero  with  the  exception  of  (XHAT(3)  -  XX(3))  ?/hich  remai- 
ned constant  at  about  .035.  Again  it  seems  that  the 
identification  of  x(3)  will  be  dificult  to  obtain. 

Next  an  experiment  was  performed  with  Eq.  (6.20) 
with  r/q=  1,  (XHATK  =   (XX) ^   but  P  -  (si   which  diver- 
ged   again  confirming  the  influence  of  the  initial  con- 
ditions when  combined  with  noise. 

The  next  set  of  experiments  were  performed  with 
the  noise  level  at  1  ,  as  before,  but  r/q  -   2,  5  and 
10. 

Comparing  the  results  from  Equations  (6.2l)  for" 
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r/q      1    and   r/q      2,    (XHAT)Q  =  [o]  and   P^   =  [2],    In   the    later 
case    the    estimate    does    not  blow   up  although   the   error 
on  the   estimate    and   the   variance    both  Increase    for  x(3) 
and  x(4),    while    for  x(5)    the   error   on  the    estimate    and 
the   variance   decreases,   which  demonstrates    some    Impro- 
vement.   The    cases  r/q  =  5  r/q   =  10  represent  no   advan- 
tage   over   r/q  =2,    presenting   the    same    trends. 

Similar    simulation  was   performed  with  Equations 
(6.20).   Abain  an  improvement  Is   noted  when  r/q   goes 
up  to    two  but   no   further   improvement  is   obtained  when 
v/q  =   5   or   10.    On  the    contrary,    therdecreasing  in  the 
error   seems   to   take    place   more    slowly   in  the    later 
case.    In  the    case   r/q  =.  2   the    advantage    of  considering 
the   non-linear   terms  when  identifying  the    linear  para- 
meters  i3   noted  because   P(3,3),    which   Increased   in  the 
case    of  Equations    (6.21),    is    seen  to  oscilate   between 
approximately  2.02   and   2.16,    and   P(4,4),    after   increa- 
sing for   the    first   couple    of  measurement,    markedly 
decreases.    The    foregoing  results  were    then  compared 

with  the    cas3    of  P   =   f5]  ,   with  clear   advantage    for    the 

— o 

former. 

The    foregoing  seems   to   indicate    that   the   values 
r/q  =    2    and   F '  =  [2j     are    the   best   suited   among  all 
other   combinations    tested.    At   this    stage    small  varia- 
tions  could   be    tried    to   further   improve    the    convergence. 

But   since    the    caracteristics   of  the   data   are    not 
expected    to    follow  very  closely    the   caracteristics   of 
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the  simulation,  it  seems  a  better  idea  to  try  to  iden- 
tify the  linear  parameters  with  the  data  obtained  from 
the  experiments  and  then  vary  the  values  of  R,  Q  and  Pq 
in  order  to  obtain  a  better  convergence.  But  before 
that, one  last  simulation  was  tried  for  the  same  values 
of  Q,  R  and  P^,  but  for  values  of   (XHAT)   closer  to  the 
true  values  :  (XHAT)  -   [0,  0,  -.2,  -.5,  .15]   .(  about 
l/2  of  the  true  values  )•  The  number  of  measurements 
was  also  increased  to  1000,  It  was  noted  that  : 

i)  P(5,5)  decreased  rapidly  to  a  value  of  .4 
(after  10  measurements  )  and  then  decreased  oscillating 
until  it  reached  a  value  of  ,057.  The  error  on  the  esti- 
mate, after  1000  measurements,  was  oscillating  between 
±.1  . 

ii)  P(4,4)  decreased  less  rapidly  to  a  value  of 
•566  (  after  30  measurements  )  and  then  jumped  to  a 
value  of  .06  where  it  kept  oscillating  for  about  350 
measurements  decreasing  finally  to  a  value  of  ,025  where 
it  stayed  oscillating  until  the  end.  The  error  on  the 
estimate  oscillated  between  ±  1  for  400  measurements 
and  then  decreased  to  zero  oscillating  about  that 
point. 

iii)  P(3,3)  increased  to  a  value  of  2.22  after 
400  measurements,  oscillated  for  a  while  and  finally 
decreased  to  reach  at  the  end  a  value  of  1.29,  with, 
tendency  for  further  decrease.  The  error  on  the  estimate 
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increased  to  a  value  of  3.3  after  420  measurements, 
oscillated  about  that  point  for  a  while  and  finally 
decreased,  oscillating  to  a  value  of  1  at  the  end. 

Note  that  the  turnover  point  seems  to  be  the 
400th  measurement  and  that  a  set  of  initial  conditions 
on  the  estimate  closer  to  the  true  value  really  impro- 
ves convergence.  This  fact  confirms  the  believe  that, 
contrary  to  the  linear  case,  the  initial  conditions  on 
both  the  estimate  and  its  covariance  play  an  important 
in  the  success  or  failure  of  the  identification. 

2 .  Identification  Using  Experimental  Data 

Following  the  conclusions  of  last  section 

a  run  was  made  with  Q  =  [j»l]  ,  R  =[.2]  and  P  —  [2.]  • 

The  set  of  initial  conditions  was  as  prescribed  by 

Equations  (6.17).  Simultaneously,  runs  with  the  values 

of  Q,  R  and  P   as  indicated  below  were  also  made.  The 
— *  —     — o 

set  of  initial  conditions,  however,  was  not  changed  at 
this  stage.  It  seemed  preferible,  since  the  above  set 
is  believed  to  be,  among  the  available  predictions,  the 
closest  to  the  true  value,  to  find  first  the  best  set 
of  Q,  R,  and  P   changing  (XHAT)   afterwards. 
Values  of  Q,  R  and  PQ  : 

R  =  [.1]  ,   [.2]  ,  [.3] 

2o  =  [1.]  ,  [2.]  ,  [5.] 
in  the  following  combinations  of  R  and  PQ  : 
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(  [.05],  [2.] 
(  [.1],[2.J 
(  [.8  ],[!.] 
(  [.3],[l.j 
For   the    group    (     [, 


,     (     f.l],[5.]    ) 

,     (      [.2],  [2.]    )     ,     (      [,2]    ,[5.J     ) 


>2  ]  »C2»]   )    the   results   obtained 
with  the    first  measurements  were   more    or   less   as 
expected.    P(5,5)    decreased  rapidly,    P(4,4)    decreased 
less   rapidly  and   P(3,3)    increased   steadily.    Then  P(4,4) 
came    close    to   zero   and   be c ame   no g a . t i ve    for   a  while,    then 
crossed   zero   again  and  went  up   to   2,48,    decreased   to 
•270   and   increased   again  to    the    final  value    of   •421. 
Also   at    the    end   P(5,5)    =   .059   and   P(3,3)  -    2,473    . 
On   the    other   runs   a   similar  behaviour  was   noted.   For 
R-  [.2]    and   P  =  [l»J    »    P(4,4)    crossed   zero   and  kept 
decreasing  up   to   a   final  value    of  -1,589.    P(5,5) 
decreased    to   reach   the   value    of    .064    (    at   the    same 
point   at  which  P(4,4)    became   negative    )    and   then 
increased    to   reach  a  value    of    .123    at   the   end,    P(3,3) 
increased   steadily   to   reach   1.462   at   the   end.   For   the 
case    R-[,2]  ,    P0=rr5.1      it  was  even  worse    since    at 
the   end   P(3,3)=    -43,225        and   XHAT(3)  =    252,70,*    For 
R-|.l|   ,    P^—  |5.|     ,    P(3,3)    increased    steadily,    P(4,4) 
decreased   to   a  value    of    ,173   and    then  increased    again 
to   reach   at   the   end    .425   with  XHAT  -  -14.513    (    ?    ). 
P(5,5)    decreased   to  a  value    of    .028   and  kept  oscilla  - 
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ting  about  that  value,  but  with  XHAT(5)  always  Increa- 
sing. For  R  -  [«l]  ,  P  =  [2.]   the  behaviour  was  similar 

but  not  so  extreme  due  to  the  smaller  value  of  P   • 

— o 

For  R-  [.05]  ,  ^-[2.]   ,  P(3,3)  always  increased 
while  P(4,4)  and  P(5,5)  decreased  first  to  low  values 
(  .60  and  .20  )  and  increased  afterwards,  never  becoming 
negative.  For  the  last  run  of  this  group  (  R  =  £«3l  , 
P_0=["l«    )>    the  type  of  behaviour  was  similar  to 
(  R=[.2]  ,  P  =  |_l. J   )  but  much  more  extreme  — —  the 
last  recorded  value  of   P(4,4)  was  -226.1  at  the  350th 
measurement  because  some  cycles  later  the  computer 
stopped  by  overflow. 

The  fact  that  P(4,4)  became  negative  posed  the 
question  of  the  correctness  of  the  computerized  scheme. 
However  the  program  after  thourough  checking  have  been 
working  successfully  with  linear  systems.  Besides  no 
such  strange  behaviour  was  found  during  the  simulation. 
There  fore  the  program  was  checked  again  in  detail 
but  nothing  could  be  found  able  to  explain  the  fore?--, 
going. 

It  was  sugested  then  that  as  P(4,4)  approaches 
zero,  the  computer  could  introduce  some  errors  thus 
making  it  negative.  The  process  would  be  self sustaining 
afterwards  *   Assuming  that  to  be  the  cause,  the  solution 
proposed  was,  whenever  P(I,I)  became  negative,  to 
substitute  it  by  it3  positive  counterpart.  Accordingly 
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the  following  statements  were  introduced  in  the  Kalman 
filter  (  KLMN  subroutine  )  : 
DO  21  I    1,  NX 

IF(  P(I,I).LT.0.0)  P(I,I)=  -P(I,I) 
21  CONTINUE 
The  effects  of  this  procedure  were  then  verified 
on  the  following  combinations  of  R  and  P   : 
(   [.2]  ,  [l.J  ) 
(   [-2]  ,  [5.]  ) 
(   [.3]  ,  [  1.2  ) 
For  the  second  case,  this  had  as  a  consequence  to 
hold  XHAT(4)  within  bounds  longer  than  before,  only  to 
diverge  more  violently  tovmrds  the  end  (  XHAT(4)=735.  ). 
The  ether  two  parameters  had  also  large  values*  For  the 
first  case,  the  influence  of  the  above  procedure  is  more 

evident,  since  originally  XHAT(4)  decreased,  crossed 
zero  at  approximately  the  100th  measurement  (  point 
where  P(4,4)  became  negative  )  and  went  on  decreasing. 
With  the  presence  of  the  above  statements,  it  behaved 
In  the  same  way  for  the  first  100  measurements,  then 
stayed  at  a  value  of  „3  until  the  300th  measurement, 
starting  a  slow  decrease  through  zero  to  reach  the 
value  of  -2.605  at  the  end  •  Acordingly  the  variance 

P(4,4),  after  touching  zero  increased  to  a  final  value 
of  .242  .  The  same  type  of  behaviour  is  also  noted  In 
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both  XHAT(3)    — diverging  with  negative   values   before, 

with  P(3,3)    increasing  and   being  constant   now  with 

P(3,3)    increasing   also,    and  XHAT(5)    — ■    going   away 

from  zero   with  positive   values  before    and   going  away 

from  zero   but  with   negative    values    and   at   a  much    slower 

rate,    with   the    introduction  of   the    above    statements* 

P(4,4)    seemed    to   be    bouncing  between  zero   and    some 

finite   positive   values,   but   since    the   values   of  P   were 

only  recorded  every  50   measurements,    no  valid   conclusion 

can  be    taken  about    that   fact.    P(3,3)    always    increased 

and   P(5,5)    remained   at   approximately   .057    • 

It  was   decided    then  to   modify   the    structure    of 

the    state   equations   to   check  if  another   form  of   the 

equations  would   converge    to  a   solution.    In  Equation 

(6.16),    for  example    the    t?/o  groups   forming   the 

coefficient   of     x(2)    -   u(l)      could   be    considered   as 

one    parameter,    being   the    new   forms   of     x(4)    and   x(5)    : 

K      -  VI. BG-  1 

10   x(4)      -        — ™— ~~  . — S 

p  xx 

(7.1) 
W.  BG  1 

10  x(5)      ~ .—-« 

~Kf>   4-   lxx 


CO' 


And   the    new   form  of   the    state   equations  would   then  be    : 

y(l)     =     .01  x(3)    x(l)  -f   x(4)  [x(2)    -   u(l)J 

+  10   x(5)    x(2) -f- .0126  x(l)    AES(    x(l)    ) 

3      (7.2) 

+   [  .0312   -f    e0167   x(5)]  Qx(2)    -   u(l)  " 

-    .0167    x(5)    x(2)5 
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y(2)  =  x(l) 

For  comparison  only  the  following  valued  of  R  and 

P  were  considered  : 
— o 

i)   R  =  [.2  ]  and  P^  [5.] 
li)  R=[.3]  and  P0=[l.l 

Also  and  to  check  the  behaviour  of  the  system  at 
high  values  of  Q  and  R,  the  following  run  was  made  : 

iiD  a=[i-j     »  R=[2«]  >   Io=[3--] 

The   main  caracteristics   of  these   runs  v/as   that   the 
parameters  did   not  oscillate,    or   oscillated  with  a  very 

long  period.    Only  in  the    case    R-[\2]      ,   Zo^C5]       » 
XHAT(3)    went   over      1,      at    the    310th  measurement, 
approximately.    XHAT(4)    and  XHAT(5)    which   started   at 
-.373   and   +.274  respectively,    crossed  each  other   and 
zero   at   the    520th  measurement  reaching  at   the   end   the 
values    +.359    and  -.422      and    showing   in  between  a 
behaviour      completely  linear.    P(3,3)    increased    to   a  value 
of  5.034   at   the    100th  measurement   and   then  decreased   as 
steadily  as    it  had    increased    to    a    final  value    of  3.412* 
Similarly,    P(4,4)    increased   to    a  value    of  5.003    in  the 
first   15  measurements   and   then  decreased   to   a   final 
value    of  4. 338,    P(5,5)    from  the    start  decreased   oscilla- 
ting until   it  reached   a   final  value    of    .384    •    This 
looked   like    being     the    object   of   the    search  apart   from 
the    fact   that   the   values    the    parameters   reached   at   the 
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end  were    not   to  be   expected   on  physical   grounds, 
specially     XHAT(5)       (    this  would  mean  that   the    center   of 
bouyancy  had    to  be    above    the    center   of  gravity  which   is 
not   true    from  the    calculated   curves    ).    A.  possibility   to 
consider   is    that    the   variables   go    through  a  cycle    of 
very  long  period   and   therefore    the    convergence    will   take 
place    later   in  the    process.    For    the    case    R      [•3~]      , 
P        [1»1         *    only   the   value    of  P(5,5)    decreased   and   the 

parameters   changed  their  value   very   slowly  and   therefore 
nothing  could   be    said   about   such  a  process.    It   could  be 
an  oscillatory  motion  of   a  period  even  longer   than 
before    or   it  could   be    a   steady,    but  very   slowly  divergent 
process.    In  the    last   case    considered,    the   parameters  did 

not   change    but   the    diagonal  elements   of   the    covariance 
matrix  increased   and  quite   rapidly  to   reach  values   of 

[1.91    ,    2.239    ,    43.289    ,    48.195    ,    32.620]      .    Note    that 
that   the   variance    of  XHAT(l)    and  XHAT(2)    was,    for    all 
the    cases    considered,    simulation  included,    of   the    order 
of   .025   and    .012   after   the    first  20    to  30  measurements. 
Therefore    this    last   case    is   of  no   consequence » 

In  summary   in   this    chapter    the    results    from   the 
identification  were    presented,    both   for   the    simulation 
carried   out   and   for  use    of   the    data.    It  was   pointed  out 
of   the    difficulty  of   the    system  to   converge.    A    se t   of 
values    of   the   different  covariance   matrices   was    selected 
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frora   the    simulation  studies    as   being   those    that  offered 
more    possibilities   fthen  using  real   data.    However    the 
results   didn't  meet   the   expectations,    being  necessary 
to  proceed   in  looking   a   set   that  would  make    the    system 

converge    •    In  stead   it  was   found    that   some    of   the 
diagonal   terms   of  P  would  become    negative.    First   it  was 
thought   to  be    resultant   from  a  mistake    in   the    compute- 
rised     scheme    and    later    to  errors    introduced  by   the 

computer    itself  and   corrected   accordingly.    Even  then 
no   solution  could   be    found   after    1000   measurements    • 
Then  another   form  of   the    state   equations   was   obtained 
by  combining   two   of   the    parameters.    This    later    system  of 

equations   presented    some    interesting  caracteristics, 
but  nothing  could  be    said   in  definite. 
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CHAPTER  8 

CONCLUSIONS   AND   RECOMENDATIONS 

From  Chapter   7   it  was   quite    clear   that   convergence 
of    the   estimate    of   the    system  to   its    true    value    was,    in 
general,    not  possible    to   obtain.    Reasons    that  could 
account  for   this   failure   in  convergence   were    then  pointed 
out   as   follows    : 

1  -    Possible   mistakes   in  the    computerized 
scheme • 

2  -   Not   accurate    prediction  of   the    initial 
conditions   of Uthe ! system. 

3  -   Data   that  measured   the    dynamical  behaviour 
was   not  accurate   enough. 

4  -   The   input   to   the  model  didn't  really  hatfe 
the    same   caracteristics   as   the   input   to   the    system 

in  the   experiments. 

To   those    reasons    above    another   may  be    pointed 
out   that  could  be    partly  responsible    for   the    non- 
convergence    of   the    system,    namely   : 

5  -   Incompatibility  between  the    system 
structure    and   the   particular    technique    that  was   chosen   • 

With  respect   to   the    first  reason,    the    fact    , 
mentioned   already   that    this    same    procedure   was   used 
before    in   the    identification  of  linear    static    systems 
successfully,    although  not   a   suficient  reason  to 
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completely  exclude  this  possibility,  at  least  decreases 
the  probability  of  such  an  event.  Therefore  it  is  recom- 
mended that  if  future  use  is  intended  to  such  a  formu- 
lation, the  program  be  rechecked  for  mistakes. 

The  not  accurate  prediction  of  the  Initial  condi- 
tions can  be  a  reason  of  some  standing  if  the  methods 
used  for  prediction  fall  under  the  category  of  empirical 

or  semiempirical  formulations  from  which  the  final 
result  have  not  the  reliability  required*  For  the  calcu- 
lation of  the  motion  parameters  for  roll,  for  example, 
the  three  categories  in  which  Chapter  3  was  divided 
also  provide  a  good  basis  for  classification  of 
accuracy,  starting  with  the  hydrostatic  coefficients 
at  the  top  and  the  non- linear  coefficients  at  the 

bottom  of  the  scale.  It  is  believed  that  for  the  case 
of  excitation  due  to  regular  waves  a  formulation  as 
presented  by  Ref.  IS  and  19  for  linear  added  mass, 
added  moment  of  inertia  and  damping  coefficient   is 
completely  adequate  in  the  sense  that  the  system  will 
converge  if  the  values  of  those  parameters  constitute 
the  only  approximation  involved.  This  cinclusion  Is 
reached  with  information  on  the  simulations  carried 
out  in  Chapter  7. 

V.Tien  trying  to  digitize  the  data  from  the  towing 
tank  experiments,  it  was  believed  that  this  item 
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played   a  major   role    in   the    success   or    failure    of  the 
identification.    This   is   not   so   anymore,    when  one    compares 
the    data   from  the    tank  "with   the    measurements   generated 
in  the    simulation,    where    the    "    level   ,!  i  of  noise    is  much 
higher   than  in  the    former   case.    This   is    not   the    same    as 
saying  that  it   is   not  necessary   to   care    about   the 
experiments.    If  nothing  else    a  careful  planning     will 
save    a   lot   of   time    and   labor,    waisted   otherv/ise    trying 
to   correct  past  mistakes* 

The    accurate    determination  of   the    input   to   the 
system  as   a   function  of  time    seems    to   play  an  important 
role    in   the    success    of   the    identification.    Intuitively 
if  the    behaviour  of   two   systems' is: to   be    compared,    it 

is  necessary  that   the    input  be    the    same.  Unfortunately 
this  was   not   the   case.    Both  methods   presented  in 

Chapters   3   and  4    (    using  the   Haskind  Relations   or 

assuming  the   major   part   of   the   excitation  to  be   due    to 
changes   in  bouyancy  due    to: the   passage    of  the  -.wave    ), 
will   not   give   very  accurate   results,    the    first     because 
it  does   not  present   the   instaneous  values   of   the    input, 
but   an  expression  that   assumes   before   hand    an  harmonic 
tims   dependency,    the    second  because,    although  it 

expresses  the  input  as  the  instantaneous  function  of 
a  quantity  being  measured  in  the  experiment,  exprimes 
a  very  rough  approximation  of   such  dependence.    However 
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the  second  was  used  in  the  identification  procedure  and 
it  is  still  believed  that,  por  this  particular  case  the 
second  yelds  better  results. 

The  incompatibility  between  the  system  structure 
and  the  method  of  identification  exists  because  some  of 
the  assumptions  set  forth  in  the  derivation  of  this 
particular  method  of  identification  were  not  met  by  the 
system.  One  of  such  assumptions   was  that  the  noise  in 
the  system,  both  environmental  and  measurement  noise 
were  uncorrelated  with  each  other  and  with  the  state 
of  the  system  (  Equations  6.5  and  6.6  ).  This  last 
point  was  not  met  by  the  vehicle  structure  since,  for 
example,  finite  water  depth  effects,  interference 
between  the  vehicle  and  the  walls,  etc.,  would  in  one 
way  or  another  be  correlated  with  motions  of  the 
vehicle.  This  fact  as  a  reason  for  failure  to  achieve 
convergence  gained  importance  when  comparing  the 
results  of  the  simulation  with  the  results  of  the 
identification,  where  the  former,  being  bad,  were  a 
bit  more  encouraging  than  the  last. 

With  the  above  in  mind,  it  seems  that  in  order  to 
improve  the  application  of  the  identification  techniques 
to  problems  of  the  type  discussed  herein,  it  will  be 
necessary  : 

1  -  Generate  a  vehicle  environment  that  can  be 
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better  approximated  as  an  uncorrelated  random  process. 
This  sugests  the  possibility  of  using  random  seas  to 
generate  the  data  for  the  identification. 

2  -  Use  another  method  of  system  identification, 
such  as  the  model  re f fere nee  technique  briefly 
described  in  Chapter  6,  which  do  not  make  such  restric- 
tive assumptions. 
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APPENDIX  1 
HYDROSTATIC  COEFFICIENTS 

1  -   DISCRIPTION  OF  THE  COMPUTER  PROGRAMS  TO  OB- 
TAIN A  POLYNOMIAL  SURFACE  REPRESENTATION 

2  -   VALUES  OF  THE  HYDROSTATIC  COEFFICIENTS  FOR  A 
MODEL  OF  THE  "  MARINER  n   HULL  FORM 
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VALUES   OF   THE   HYDROSTATIC    COEFFICIENTS 
Since    the    identification   of   the    motion  parameters   in 
coupled  heave    and  roll  was   not  performed,    it  was   necessary 
to   calculate    the   hydrostatic   coupling   coefficients.    To   cal- 
culate   only   the    coefficients   Zz    ,    Zz2,    K^    and  103   a  much 
simpler  method    than  resorting   to    the   Kerwin  polynomial  re- 
presentation  and    that  will   probabilly  ye Id  more    accurate 
results   can  be   used    •    The   reason  for   such  is    that  in  com- 
puterized  solutions  using  least   squares,    it   is   not  usual- 
ly possible    to   specify   some    a  priori   knowledge    that   one 
might  have    about   the    surface    or    the    curve.    For  example    in 


roll  K(^=0)      and  K^2 
symmetry.    Also  K^ 


,        are    zero   due    to  port   and   starboard": 
v  0=    -  W---GM    •    A   least    square    method   gave 
nonzero  values   for   the    first    two   coef f iciente s. 

Looking   just   for    the    secon  non  zero   term  in   the    Tay- 
lor's expansion  and  knowing  the    first   it's  much  easier   to 
calculate   manually   the   missing   term.   Using  this   procedure 
and  being,    for   roll  KQ      -  W*GM      -2.159   LB»FT,    the    value    that 
will   approximate    the    moment   curve   better  will  be    : 

Between    JZ^i.175  l/3'-*KQ3      =        6.814   LB*FT 

^~ ±,350  *     -1.254    LB*FT 

<£=.+  .  525  ■=     -3.007    LB*FT 

Since    the    roll   angle   will   go   over    .175   but  does   not 
reach    .350   In  the   experiments    the   value    choose n  was: 
1/3' K03      -      6.814    LB-::-FT 

Similarly   for   heave    the    values    found  were 
Z-fc  =  -  j°*g*Alw    =     -194.709    LE/FT 

l/2'Zz2  =       -84.528    LB/FT2 
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APPENDIX  2 

LINEAR  HYDRODYNAMIC  COEFFICIENTS 

1  -   DISCRIPTION  OF  THE  COMPUTER  PROGRAMS  TO  CAL- 
CULATE ADDED  MASS,  ADDED  MOMENT  OF  INERTIA  AND  DAMPING 
COEFFICIENTS  FOR  TViO  DIMENSIONAL  CYLINDERS  IN  HEAVE,  SWAY 
AND  ROLL 

2  -   VALUES  OF  ADDED  MASS,  ADDED  MOMENT  OF  INERTIA 
AND  DAMPING  COEFFICIENTS  FOR  A  MODEL  OF  THE  "  MARINER  " 
HULL  FORM  IN  HEAVE  AND  ROLL 
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PROGRAM    INPUT 
Cards    1-4  FORMAT (20A4) 

The    purpose    of   these    cards   Is    to   provide    titles    for 
the    computed    output.   Vfhen  titles    are    not  desired  blank 
cards    should  be    provided. 
Card    5   -NOK,MAXB,NH  FORMAT (31 6) 

NOK  -   Total   number   of  frequencies,    at  which   the    cal- 
culations will  be    performed. (MAX      50) 

MAXB  -Integer   correspondent   to   the   horizontal   plane 
of  maximum  beam.    Start   the   count  at   the   keel. 

NH  -   Number   of  depths   at  which   the    calculations   will 
be    performed. (MAX      10) 
Card    6   -    CAY(K),K   1,N0K  FORMAT (5F12. 7) 

CAY  -   Nondimensional   frequencies        -  a/g 

-  radian   frequency 

a  -  half  the  beam  for  a  floating  cylinder 

-  half  the  draft  for  a  submerged  cylinder 
g  -  gravitational  acceleration 

Card  7  -(DEP(J),J  1,NH),CR        FORMAT (5F12. 7) 

DEP   -   Values   of   the    depths   of   submergence,   measured 
frogi  the    surface    to   the    top  of  the   cylinder.   For   a   floa- 
ting cylinder,    NH  1   and   DEP   0    . 
Card   8   -    MD,N0N  FORMAT (21 6) 

MD   -   Control   integer   -    1   for    a   submerged  body 

2    for   a   floating  body 
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NON  -   Number   of  segments    In  which   the    right  hand    side 
half  contour   is   divided    (MAX      45) 

CR  -  Positive  vertical  distance  from  the  center  of 
rotation  to  the  free  surface,  for  a  floating  cylinder  or 
from  the    center   of  rotation  to   the    intersection  of   the 

section  contour   with   the   plane    of   symmetry(    on   the    side 
closer   to   the    surface. 
Card    9    -   X(J),J   1,    (NON  1)  FORMAT (5F12. 7) 

X(J)    -   Horizontal  offsets   of  the   end  points   of  the 
line    segments,   measured   from  the   vertical  planr   of  symme- 
try to    the    right  hand    side    of  the    contour,    starting  at   the 
keel.    The    last  point  is    the    intersection  of   the    section 
contour  either   with  the    free    surface    (    floating  cylinder    ) 
or  with  the   vertical  plans?   of   symmetry(    submerged  cylinder    ) 
Card   10   -   Y(J),J   1,    (NON  1)  FORMAT (5F12. 7) 

Y(J)    -   Negative   vertical   offsets   of  the   end  points   of 
the    line    segments,    measured   from  th   free    surface,    for   a 
floating  cylinder    or   from  the    intersection  of  the    contour 
with    the   vertical  plane    of   symmetry,    for    a   submerged   cy- 
linder,   starting  from  the    keel*   Again  only   the    right  hand 
side    of   the    contour  will  be    considered* 
Card   11   -    MODE  F0RMAT(I6) 

MODE   -   Control   integer   -    1      -  heave 

2  -    sway 

3  -   roll 

Card  12  -  FORMAT (20A4) 
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Provides  a  title  for  the  output  of  the  program.  If 
the  title  is  not  desired,  a  blank  card  should  be  provided, 
If  more  than  one  mode  is  desired,  adicional  pairs  of 
cards  11  and  12  should  be  provided,  accordingly. 
Card  13  -  MODE  FORMAT (16) 

MODE   -   Control  variable 
0  -  End  of  data 

Negative  integer   -   A  new  geometry  follows 
starting  Card  5. 
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APPENDIX  3 

NON-LINEAR  MOTION  PARAMETERS 

1   -   DETERMINATION  OF  HYDRODYNAMIC  COEFFICIENTS 
USING  EMPIRICAL  METHODS 
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1      -      ADDED   MASS   AND  ADDED  MOMENT   OF    INERTIA 
Using  Ref.    20  p.    61   to    63        -      Heave 


Am  =  t--"-J°*   V2/       ^'7    (dx) 


(A3.1) 


where     P    -  water  density  1.94  LB---SEC2/FT4 

L  -  lenght  of  the  model       5.5  FT 
£  -  coefficient  from  fig.  17  ^   1. 
y  -  contour  of  the  waterline,  aproximated  by  a 
parabola  of  the  form 


y=   b/2 


1     172 


n 


(A3. 2) 


79   FT 


B  —   beam  of  the   model 

n    -  cv/(i-ew) 

C      -  waterplane   coefficient   -   Av/-j/B»-L      =       .724 
f     -   coefficient,    calculated  according  to 


^  =    1/    [  1  +.8*B/l] 


.395  (A3. 3) 


the  following  form  may  be  obtained 


2  .2 


A  m  =  .772*J*  p  *L*B**cy(H-Cw)    (A3. 4) 


Plugging  in  the   values 

Am  =  1.41   LB::-S3C2/FT 
Being  the   mass    of   the   vehicle      -      m     1.43   LB---SEC  /FT 
this   result  may  be    checked  by  measuring  the    natural  period 
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of  heaving  motions  and  using  results  from  the  hydrostatic 


calculations  Z_  -  -  -"-g-^A 


irl 


OJ 


n 


-    Z. 


m  +a  m 


-194.05    LB/FT 


8.29    SEC(-l) 


The   value    found   for    the    natural  undamped    frequency  was 
8.038    (SSC-1). 

Similarly  for  roll,    the   added  moment   of  inertia  is 

expressed   as   a  percentage   of  the   moment   of  inertia.    This 
is   found,    from  data  reffering  to   full   scale   vehicles,    acor« 

din/7   to   the    formula 


XXX  *     w/g*((B2*Cw)/ll.4*-Cb  -*"  D2/12    )    T0N*3EC' 
.02855    LBSEC2 
From  fig.      24   of  the    above   re f fere nee 
aIxx=   .286-"-Ixx 


Therefore 


I        +     a  I 
xx  xx 


.03665   LB-"- SEC 


(A3. 5) 


From  the  hydrostatic  calculations 

K0--  W-"-GM    -     -2.159   LB*FT 
The    resulting  natural   frequency   is    then 
(Jn    -       7.7      SEG(-l) 


while    the   value    found   in   the   measurements  was 


uJ^  -  4.14  SEC(-l) 
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The  values  of  added  mass  and  of  adde  s  moment  of 
inertia,  calculated  in  Appendix  2  were  : 
A (3, 5)  =   1.24  LB-::-SEC2/FT 
A (4, 4)   =   .0176  LB-::-SECS*FT 
The  value  of  roll  added  moment  of  inertia  is  seen  to 
be  quite  different  the  results  of  both  methods  are  compared, 

2  -  DAMPING  COEFFICIENT 
For  heave  the  following  expression,  taken  from  page 
153  of  the  above  re f fere nee  was  used  : 

2 

2  y'ML.B.  C_J 

where    the   parameters   are    as  before    and   also    : 
^    -   wave    lenght   =     3.11  FT   for  u;  =  8.04   SEC(-l) 


Z^,        -1/2P2  '/ML.B.C^^TAl.iMA)  (A3. 6) 


YC  -   coefficient   taken  from    (20)    fig.    52     =      ,460 

T      -   draft  of   the   vehicle      =    .28  FT 

F£i-   ±   coefficient   taken  from   (20)fig.    53         -    .200 

For   the   values   above    :  Z  =  -2.336   LB-::-SEC/FT 

w 

while    when  calculated  by   the    methods   of  Appendix   2 

Z^  =  -7.4543   LB-::-SEC/FT  which   is    three    times   the 

value   presented   above. 

For  roll,    the   calculation  of  linear   and   non-linear 
damping  coefficients  were   carried  simultaneously,    for  each 
transverse    section  of  the   vehicle,    using  Equations   21,    22, 
23   or   24   of  page    148   of   the    above    re f fere nee,    deppending  on 
the    shape    of   the    section^    The   value    of   the    damping  for   the 
whole   vehicle   was   obtained  by   summing  the    sectional  results 
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over      the    lenght   of   the   vehicle    in  a   "    strip   theory   "   way. 

Therefore,    using  the    following  equations    the    damping  is 
calculated. 


KDi    -l/2jobiAli 


°i<»o.y»i>l*   (bi/g)frv 


p  (A3. 7) 


kdi  -*rbiAii 


Coi<VVVlpl+ 


72  (Bi/Ti-5. 07^11 


KDi    "*fbiAli 


Coi(VW|PK 


(bi/g^cJ-4    10 


36(Bi/Ti-2.) 


(A3. 8) 


(A3. 9) 


»»!  -^fV1! 


57|p|        + 


g(B1/T1-2.82")    4 


[bi/gTTlO3" 


(A3. 10) 


ft 


V 

>° 

— , 

I 


Fig.    A3.1 


where   Eq.    A3. 7    and  A3. 8   are   used   for   rectangular    sections 
of  different   value    of   r0      ri/hi    ,   Eq.    A3. 9   is   used  with 
elliptical    sections   and  Eq  A3. 10  with   terminal    section  of 
shape    close    to   triangular,    the    following  values  were 
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found   for   linear   and   non  linear   damping  coefficients: 

Kp      -       -,00786   LB*FT-» SEC  -    (x>  --  4.145.  SEC  (-1) 

=     -.06835   LB*FT* SEC  -     cj  =  8.290. SEC(-l) 

Kn      -    -.01605    LB-::- FT-»- SEC2  -      all    frequencies 

The   values    found   in  Apenndix  3    for   the    linear   damping 
coefficients   were    respectively     -.0022   LB*FT»SEC' 

-.0363    LB*FT*SEC 
For    the   nonlinear  heave   damping  coefficient,    the  va- 
lues  for  each   section  were    obtained  by  considering  each   se- 
ction    substituted  by   an  equivalent  cylinder   of  the    same 
diameter  when  the   model  moves   away   from  the   water    (    assu- 
ming the    only  drag  present   to  be    separation  drag   )    and   sub- 
stituted by   a  cylinder   of  the    same    lenght   of  wetted  con- 
tour when  the   model  moves   Into   the   water    (    assuming   the   on- 
ly drag  present   to   be    skin  friction  drag   )    •      Data   from 
infinite    cylinders  was    then  used,    the    correspondent     values 
of  drag  coefficient   averaged    over   a  cycle,    and    the    results 
summed  over    the    lenght   of  the   model.    The   velocity  entered 
in  the    calculation  of   the    Reynolds    Number  was   a  RIv.S   type 
of  velocity   taken  from  the   records   of   the   experiments. 
The    final  value    obtained  was: 


Zw|w|    =      -1.767   LB*SECC/FT' 
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APPENDIX  4 

CIRCUITRY  AND  DIGITIZING  PROGRAM 

1   -   PROGRAM  TO  DIGITIZE  DATA  AND  CORRESPONDENT 
PATCHING  IN  THE  ANALOG  COMPUTER 
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APPENDIX  5 
IDENTIFICATION  PROCEDURE 

1  -   DISCRIPTION  OP  THE  COMPUTER  PROGRAMS  USED  IN 
THE  IDENTIFICATION  PROCEDURE 

2  -   EXAMPLE  OF  THE  BEHAVIOUR  OF  VARIABLES  AND  RE- 
LATED  VARIANCES  DURING  IDENTIFICATION 
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